arXiv:astro-ph/0211649vl 29 Nov 2002 


Astronomy &: Astrophysics manuscript no. sznt'bab February 2, 2008 

(DOI: will be inserted by hand later) 


The Non-Thermal Sunyaev-Zel’clovich Effect in Clusters of 

Galaxies 

S. Colafrancesco^, P. Marchegiani^ and E. Palladino^ 

^ INAF - Osservatorio Astronomico di Roma via Frascati 33, 1-00040 Monteporzio, Italy. 

Email: cola@coma.mporzio.astro.it 

^ Dipartimento di Fisica, Universita di Roma “La Sapienza”, Piazzale A. Moro 2, I-OOO Roma, Italy 
Received 21 May 2002 / Accepted 13 August 2002 

Abstract. In this paper we provide an general derivation of the non-thermal Sunyaev-Zel’dovich (SZ) effect in 
galaxy clusters which is exact in the Thomson limit to any approximation order in the optical depth r. The general 
approach we use allows also to obtain an exact derivation of the thermal SZ effect in a self-consistent framework. 
Such a general derivation is obtained using the full relativistic formalism and overcoming the limitations of the 
Kompaneets and of the single scattering approximations. We compare our exact results with those obtained at 
different approximation orders in r and we give estimates of the precision fit. We verified that the third order 
approximation yields a quite good description of the spectral distortion induced by the Comptonization of CMB 
photons in the cluster atmosphere. In our general derivation, we show that the spectral shape of the thermal 
and non-thermal SZ effect depends not only from the frequency but also from the cluster parameters, like the 
electron pressure and the optical depth and from the energy spectrum of the electron population. We also show 
that the spatial distribution of the thermal and non-thermal SZ effect in clusters depends on a combination of 
the cluster parameters and on the spectral features of the effect. To have a consistent description of the SZ effect 
in clusters containing non-thermal phenomena, we also evaluate in a consistent way - for the first time - the 
total SZ effect produced by a combination of thermal and non-thermal electron population residing in the same 
environment, like is the case in radio-halo clusters. In this context we show that the location of the zero of the total 
SZ effect increases non-linearly with increasing values of the pressure ratio between the non-thermal and thermal 
electron populations, and its determination provides a unique way to determine the pressure of the relativistic 
particles residing in the cluster atmosphere. We discuss, in details, both the spectral and the spatial features of 
the total (thermal plus non-thermal) SZ effect and we provide specihc predictions for a well studied radio-halo 
cluster like A2163. Our general derivation allows also to discuss the overall SZ effect produced by a combination 
of different thermal populations residing in the cluster atmosphere. Such a general derivation of the SZ effect 
allows to consider also the CMB Comptonization induced by several electron populations. In this context, we 
discuss how the combined observations of the thermal and non-thermal SZ effect and of the other non-thermal 
emission features occurring in clusters (radio-halo, hard X-ray and EUV excesses) provide relevant constraints on 
the spectrum of the relativistic electron population and, in turn, on the presence and on the origin of non-thermal 
phenomena in galaxy clusters. We finally discuss how SZ experiments with high sensitivity and wide spectral 
coverage, beyond the coming PLANCK satellite, can definitely probe the presence of a non-thermal SZ effect in 
galaxy clusters and disentangle this source of bias from the cosmologically relevant thermal SZ effect. 

Key words. Cosmology: theory - Galaxies: clusters: general - Intergalactic medium - Radiation mechanism: 
thermal - Radiation mechanism: non-thermal - Cosmic microwave background 


1. Introduction 

Compton scattering of the Cosmic Microwave Background (CMB) radiation by hot Intra Cluster (hereafter IC) elec¬ 
trons - the Sunyaev Zel’dovich effect (Zel’dovich and Sunyaev 1969, Sunyaev & Zel’dovich 1972, 1980) - is an im¬ 
portant process whose spectral imprint on the CMB can be used as a powerful astrophysical and cosmological probe 
(see Birkinshaw 1999 for a review). Such a scattering produces a systematic shift of the CMB photons from the 
Rayleigh-Jeans (RJ) to the Wien side of the spectrum. 
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An approximate description of the scattering of an isotropic Planckian radiation field by a non-relativistic 
Maxwellian electron population can be obtained by means of the solution of the Kompaneets (1957) equation (see, 
e.g., Sunyaev & Zel’dovich 1980). The resulting change in the spectral intensity, Alth, due to the scattering of CMB 
photons by a thermal electron distribution can be written as 

A/tft = 2^^^^^^ythg{x) , (1) 


where x = hv/ksT^ is the a-dimensional frequency and the spectral shape of the effect is contained in the function 


9{x) 


(e^ - 1)2 f - 1 


( 2 ) 


which is zero at the frequency xq = 3.83 (or v = 217 GHz for a value of the CMB temperature Tq = 2.726 K), negative 
at X < xo (in the RJ side) and positive at x > xq (in the Wien side). The Comptonization parameter, yth, due to the 
thermal SZ effect is given by 


yth = [ M rieksTe , (3) 

rrieC^ J 

where Ue and Te are the electron density and temperature of the IC gas, respectively, ctt is the Thomson cross 
section, valid in the limit ^ Tq, ks is the Boltzmann constant and WeC^ is the rest mass energy of the electron. The 
Comptonization parameter is proportional to the integral along the line of sight (. of the kinetic pressure, Pth = riefcsTe, 
of the IC gas. Thus, the previous Eq.(^) can be written as 

yth = f d£ Pth , (4) 

rrieC^ J 

where the relevant dependence from the total kinetic pressure, Ptt, of the clusters along the line of sight i appears. 
The previous description of the thermal SZ effect is obtained under the Kompaneets (1957) approximation and in the 
single scattering regime of the true photon redistribution function. As such, it only provides an approximation of the 
SZ effect in galaxy clusters for low temperatures (Te ^ 3 keV) and low optical depth (r = ut / dlue ~ 10“^). 

However, recent X-ray observations have revealed the existence of many high-temperature clusters (David et al. 
1993; Arnaud et al. 1994; Markevitch et al. 1994; Mushotzky and Scharf 1997) with ksTg up to ^ 17 keV (see, e.g.. 
Tucker et al. 1998). Thus, the calculation of the thermal SZ effect from these hot clusters requires to take into account 
the appropriate relativistic corrections (see Birkinshaw 1999 and references therein). 

Analytical expressions for the SZ effect in the relativistic case have been considered by various authors (Stebbins 1997; 
Itoh et al. 1998; Challinor & Lasenby 1998; see also Birkinshaw 1999 and references therein). Some of these calculations 
(out of the Monte Carlo simulations) are still approximate since they have been carried out in the following limits: i) 
single scattering of CMB photons against the IC gas electrons; ii) diffusion limit in which the use of the Kompaneets 
equation is justified. Itoh et al. (1998) obtained higher-order relativistic corrections to the thermal SZ effect in the form 
of a Fokker-Planck expansion (generalised Kompaneets equation) and by direct integration of the Boltzmann collision 
term. Such an integration has been carried out in the full relativistic regime. So, the results of Itoh et al. (1998) can be 
regarded as exact in the framework of the single scattering approximation. Analytical fitting formulae of such derivation 
can be found in Nozawa et al. (2000) and Itoh et al. (2002). Such a formalism seems to offer a detailed description 
of the thermal SZ effect for keTe ^ 15 keV, while Monte Carlo simulations describe more correctly the thermal SZ 
effect even for fc^Tg ^ 20 keV (see, e.g., Challinor and Lasenby 1998). Recently, Dolgov et al. (2001) proposed an 
approach based on an analytical reduction of the collision integral which is contained in the Boltzmann-like collision 
equation for the SZ effect produced by a thermal population of electrons. They found a numerical solution for the 
thermal SZ effect which is valid for generic values of t and Tg. Sazonov and Sunyaev (2000) presented a derivation 
of the monochromatic redistribution function in the mildly relativistic limit which considers also quantum effects and 
the use of the Klein-Nishina cross-section. However, they still consider only the single Compton scattering limit and 
the relativistic corrections up to some intermediate order due to low-energy photons and relativistic electrons. Such 
derivation of the monochromatic redistribution function in the limit hv <C ksTe reproduces the results of Fargion et 
al. (1997). More recently, Itoh et al. (2001) also presented a calculation of the thermal SZ effect which considers the 
contribution from multiple scattering in the relativistic limit. These last authors concluded that the multiple scattering 
contribution to the thermal SZ effect is negligible in galaxy clusters. 

Another general assumption which is made in the calculation of the SZ effect is the use of a single population 
of thermal electrons which constitutes the hot (with temperature Tg ^ 10^ — 10® K), optically thin (with density 
rig ^ 10“® — 10“^ cm“®, and size R ^ a few Mpc) Intra Cluster Medium (hereafter ICM). This assumption is based 
on the evidence that the IC gas is mainly constituted by thermal electrons (and protons) which are responsible for the 
X-ray emission observed in many clusters through thermal bremsstrahlung radiation (see Sarazin 1988 for a review). 
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Nonetheless, in addition to the thermal IC gas, many galaxy clusters contain a population of relativistic electrons 
which produce a diffuse radio emission (radio halos and/or relics) via synchrotron radiation in a magnetized ICM (see, 
e.g., Feretti 2000 for a recent observational review). The electrons which are responsible for the radio halo emission 
must have energies ^ a. few GeV to radiate at frequencies ^ 30 MHz in order to reproduce the main properties 
of the observed radio halos (see, e.g., Blasi & Colafrancesco 1999; Colafrancesco & Mele 2001, and references therein). 
A few nearby clusters also show the presence of an EUV/soft X-ray excess (Lieu et al. 1999, Kaastra et al. 1999, 
2002; Bowyer 2000) and of an hard X-ray excess (Fusco-Femiano et al. 1999-2000; Rephaeli et al. 1999; Kaastra et 
al. 1999; Henriksen 1999) over the thermal bremsstrahlung radiation. These emission excesses over the thermal X-ray 
emission may be produced either by Inverse Compton Scattering (hereafter ICS) of CMB photons off an additional 
population of relativistic electrons or by a combination of thermal (reproducing the EUV excess. Lieu et al. 2000) 
and suprathermal (reproducing the hard X-ray excess by non-thermal bremsstrahlung, Blasi et al. 2000, Dogiel 2000, 
Sarazin & Kempner 2000) populations of distinct origins. However, since the inefficient non-thermal bremsstrahlung 
mechanism would require a large energy input and thus imply an excessive heating of the IC gas - which is not observed 
- a more complex electron population is required to fit both the radio halo and the EUV/hard X-ray spectra of galaxy 
clusters (Petrosian 2001). The high-energy part {E ^ 1 GeV) of such spectrum contains non-thermal electrons and 
produce synchrotron radio emission which can fit the observed radio-halo features: this is the best studied region of 
the non-thermal spectrum in galaxy clusters. 

Indeed, in many of the clusters in which the SZ effect has been detected there is also evidence for radio halo sources 
(see Colafrancesco 2002). So it is of interest to assess whether the detected SZ effect is, in fact, mainly produced by 
the thermal electron population or there is a relevant contribution from other non-thermal electron populations. We 
will show specifically in this paper that each one of the electron populations which reside in the cluster atmosphere 
produces a distinct SZ effect with peculiar spectral and spatial features. 

The description of the non-thermal SZ effect produced by a single electron population with a non-thermal spectrum 
has been attempted by various authors (McKinnon et al. 1991; Birkinshaw 1999; Ensslin & Kaiser 2000; Blasi et al. 
2000). Several limits to the non-thermal SZ effect are available in the literature (see, e.g., Birkinshaw 1999 for a review) 
from observations of galaxy clusters which contain powerful radio halo sources (such as A2163) or radio galaxies (such 
as A426), but only a few detailed analysis of the results (in terms of putting limits to the non-thermal SZ effect) have 
been possible so far. Only a single devoted search for the SZ effect expected from a relativistic population of electrons 
in the lobes of bright radio galaxies has been attempted to date (McKinnon et al. 1991). No signals were seen and a 
detailed spectral fit of the data to separate residual synchrotron and SZ effect signals was not done, so that the limits 
on the SZ effect do not strongly constrain the electron populations in the radio lobes. Also the problem of detecting the 
non-thermal SZ effect in radio-halo clusters is likely to be severe because of the associated synchrotron radio emission. 
In fact, at low radio frequencies, such a synchrotron emission could easily dominate over the small negative signal 
produced by the SZ effect. At higher frequencies there is in principle more chance to detect the non-thermal SZ effect, 
but even here there are likely to be difficulties in separating the SZ effect from the flat-spectrum component of the 
synchrotron emission (Birkinshaw 1999). 

From the theoretical point of view, preliminary calculations (Birkinshaw 1999, Ensslin & Kaiser 2000, Blasi et 
al. 2000) of the non-thermal SZ effect have been carried out in the diffusion approximation (r <C 1), in the limit 
of single scattering and for a single non-thermal population of electrons. Specifically, Ensslin & Kaiser (2000) and 
Blasi et al. (2000) considered the SZ effect produced, under the previous approximations, by a supra-thermal tail 
of the Maxwellian electron distribution claimed to exist in the Coma cluster and concluded that the SZ effect, even 
though of small amplitude, could be measurable in the sub-mm region by the next coming PLANCK experiment. 
However, Petrosian (2001) showed that the suprathermal electron distribution faces with several crucial problems, 
the main being the large heating that such electrons would induce through Coulomb collisions in the ICM. The large 
energy input of the suprathermal distribution in the ICM of Coma would heat the IC gas up to unreasonably high 
temperatures, ksTe ~ 10^® K, which are not observed. In addition, Colafrancesco (2002) noticed that dust obscuration 
does not allow any detection of the SZ signal from Coma at frequencies ^ 600 GHz. 

Matters are significantly more complicated if the full relativistic formalism is used. However, this is necessary, since 
many galaxy clusters show extended radio halos and the electrons which produce the diffuse synchrotron radio emission 
are certainly highly relativistic so that the use of the Kompaneets approximation is invalid. Moreover, the presence of 
thermal and non-thermal electrons in the same location of the ICM renders the single scattering approximation and 
the single population approach unreasonable, so that the treatment of multiple scattering among different electronic 
populations coexisting in the same cluster atmosphere is necessary to describe correctly the overall SZ effect. 

Here we derive the spectral and spatial features of the SZ effect using an exact derivation of the spectral distortion 
induced by a combination of a thermal and a non-thermal population of electrons which are present at the same time 
in the ICM. Lately, we also consider the case of the combination of several non-thermal and thermal populations. 
The plan of the paper is the following: in Sect. 2 we will provide a general derivation of the SZ effect in an exact 
formalism within the framework of our approach, i.e., considering the fully relativistic approach outlined in Birkinshaw 
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(1999) and the effect induced by multiple scattering. We work here in the Thomson limit hiy <C ruef?- Within such 
framework, we derive the exact SZ effect for a single electron population both in the thermal and non-thermal cases. 
In Sect. 3 we derive the SZ effect produced by a single population of non-thermal electrons providing both exact 
and approximate (up to any order in r) expressions for the non-thermal SZ effect. In Sect. 4 we derive the total SZ 
effect produced by a combination of two populations: a distribution of thermal electrons - like that responsible for the 
X-ray emission of galaxy clusters - and a non-thermal electron distribution - like the one responsible for the radio halo 
emission which is present in many galaxy clusters. Our general approach allows to derive the expression for the SZ 
effect produced by the combination of any electronic population. So, in Sect. 5 we consider also the SZ effect generated 
by the combination of two different thermal populations. In Sect. 6 we discuss the spatial features associated to the 
presence of a non-thermal SZ effect superposed to the thermal SZ effect. In Sect. 7 we derive limits on the presence 
and amplitude of the non-thermal SZ effect discussing specifically the cases of a few clusters in which there is evidence 
for the presence of non-thermal, high-energy electrons. We show how the possible detection of a non-thermal SZ effect 
can set relevant constraints on the relative electron population. We summarize our results and discuss our conclusions 
in the final Sect. 8 . We use Hq = 50 km s“^ Mpc“^ and flo = 1 throughout the paper unless otherwise specified. 


2. The SZ effect for galaxy clusters: a generalized approach 

In this section we derive a generalized expression for the SZ effect which is valid in the Thomson limit for a generic 
electron population in the relativistic limit and includes also the effects of multiple scattering. First we consider an 
expansion in series for the distorted spectrum I(x) in terms of the optical depth, r, of the electron population. Then, 
we consider an exact derivation of the spectral distortion using the Fourier Transform method already outlined in 
Birkinshaw (1999). 

An electron with momentum p = Pj, with P = vfc and 7 = Eejmef? increases the frequency v oi a. scattered 
CMB photon on average by the factor t = v' jv = | 7 ^ — i, where v' and v are the photon frequencies after and before 
the scattering, respectively. Thus, in the Compton scattering against relativistic electrons (7 ^ 1) a CMB photon is 
effectively removed from the CMB spectrum and is found at much higher frequencies. We work here in the Thomson 
limit, (in the electron’s rest frame ^hv <C rrieC^), which is valid for the interesting range of frequencies at which SZ 
observations are feasible. 

The redistribution function of the CMB photons scattered once by the IC electrons writes in the relativistic limit as. 


Pi{s) = / dpfe{p)Ps{s;p) 


( 5 ) 


where fe{p) is the electron momentum distribution and Ps{s;p) is the redistribution function for a mono-energetic 
electron distribution, with s = ln{t). An analytical expression for the redistribution function of CMB photons which 
suffer a single scattering, Ps{s;p), has been given by Ensslin & Kaiser (2000). The expression for such a redistribution 
function has been also given analytically by Fargion et al. (1997) and by Sazonov & Sunyaev (2000). Once the function 
Pi{s) is known, it is possible to evaluate the probability that a frequency change s is produced by a number n of 
repeated, multiple scattering. This is given by the repeated convolution 

/ -I-00 

dsi. . . dSn-lPl{si) . . . Pi{Sn-l)Pl{s - Si - ... - S„_i) = Fi(s) O ... O Pi{s) (6) 

-00 V ^ 

n times 

(see Birkinshaw 1999), where the symbol O indicates each convolution product. As a result, the location of the 
maximum of the function Pn{s) moves towards higher values of s for higher values of n and the distribution Pnis) 
widens a-symmetrically towards high values of s giving thus higher probabilities to have large frequency shifts. The 
resulting total redistribution function P{s) can be written as the sum of all the functions Pn{s), each one weighted by 
the probability that a CMB photon can suffer n scatterings, which is assumed to be Poissonian with expected value r: 


Pis) = 


-1-00 __ , 

n=0 


nl 


■Pnis) 


Pois) + TPiis) + -T^P2(s) + ■ • ■ 

5(s) -f rPiis) + ^r'^Piis) O Pi(s) 

The spectrum of the Comptonized radiation is then given by 

(•+00 


/ +00 

ds /o(a;e“®)P(s) 

-OC 


( 7 ) 


( 8 ) 
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where 


/o(x) =2 


{keTof 


{hcY — 1 


( 9 ) 


is the incident CMB spectrum in terms of the a-dimensional frequency x. 

In the following we will derive the expression for the distorted spectrum using first an expansion in series of r and 
then the exact formulas obtained with the Fourier Transform (FT) method. 


2.1. High order T-expansion 

In the calculation of the SZ effect in galaxy clusters it is usual to use the expression of P{s) which is derived in the 
single scattering approximation and in the diffusion limit, r <C 1. In these limits, the distorted spectrum writes, in our 
formalism, as: 


I{x) = Jo{x) + T 
where 


Ji(x) - Jo{x) 


/ -hOO 1*^00 

Io{xe~‘‘)PQ{s)ds = / Io{xe~'')S{s)ds = Iq{x) , 
-oo J — oo 

/ + 00 

Io{xe~^)Pi{s)ds . 

-oo 


( 10 ) 

( 11 ) 

( 12 ) 


To evaluate the SZ distorted spectrum I{x) up to higher order in r, we make use of the general expression of the series 
expansion of the function P{s) 


+ 00 


Pis) = 


n—0 


which can be written, using Eq.(^, as 

+°° ( ^\k 

k'—O 


(13) 


(14) 


The general n-th order term is obtained by selecting the terms in the double summation which contain the optical 
depth r up to the n-th power. These terms are obtained for k' = n — k, and provide the following expression for the 
series expansion coefficients: 


/_ -t\k 1 






(-l)'=P„_fe(s). 


(15) 


Inserting the coefficients a„(s) given in Eq.(15) in Eq. (H) and this last one in Eq.(8), the resulting spectral distortion 
due to the SZ effect can be written as: 


+ 00 


Hx) = ^n(a;)r” , 

n—0 

where the coefficients bnix) are given by 

k=0 ^ ' 

and 

/ + 00 

Iq(^XG (5)c?S. 

-oo 


(16) 


(17) 


(18) 


In this general approach, we can write the corrections to the distorted spectrum I{x) up to any order n in the quantity 

T. 
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2.2. Exact derivation of the SZ effect 


The redistribution function P{s) can also be obtained in an exact form considering all the terms of the series expansion 
in Eg. (|l^ . In fact, since the Fourier transform (hereafter FT) of a convolution product of two functions is equal to 
the product of the Fourier transforms of the two functions, the FT of P{s) writes as 


P{k) = e" 
where 


l + rP,{k) + K^P^{k) + ... 


— g—rgTPi(fc) _ g-T[l-Pi(fe)] 


^+oo 

Pl(*) = / Pl{s)t 


— iki 


ds 


(19) 


( 20 ) 


(Taylor & Wright 1989). The exact form of the Comptonized spectrum I(x) is then given by Eq.(g) in terms of the 
exact redistribution function 


p(,) = 






P{k)P'^^dk 


( 21 ) 


which is obtained as the anti Fourier transform of P{k) given in Eq.(|lS|). To compare the exact calculations of the SZ 
spectral distortion with those obtained in the non-relativistic limit [see eqs.(l- 2 ) as for the thermal case], it is useful 
to write the distorted spectrum in the form 


AT/ \ r,ikBTo)^ ^ 

^ (feg )2 > 

where AI{x) = I{x) — Io{x). The spectral shape of the SZ effect is contained in the function 
' AI\ 1 x^ Ai{x) 


Kx) = 


y — 1 
jhc) 


y 


( 22 ) 


(23) 


where Ai = The Comptonization parameter y is defined, in our general approach, in terms of the pressure 

P of the considered electron population: 


y = 


Ct 

rrieC^ 


Pd£ . 


(24) 


2.3. The case of a single thermal electron population 

For a thermal electron population in the non-relativistic limit with momentum distribution fe,th cx p^exp{—r]yA+]P) 
with p = rrieC^/ksTe, one writes the pressure as 


Pfh — riekBTe 

and it is easy from Eq.(p^ to re-obtain the Compton parameter in Eq.(3) as 


Uth — 


(7t 

rrieC^ 


neksTedt = T 


kaTf., 


rrieC^ 


(25) 


(26) 


(we consider here, for simplicity, an isothermal cluster). The relativistically correct expression of the function g{x) for 
a thermal population of electrons writes, at first order in r, as: 


aix) = 


_ T[h - jo] _ nieC^ 


yth 


.ksT^ 


ksTe 


b'l - jo ], 


(27) 


meC-^ 

/ 7 \ 2 

where ji = Ji functions Ji are given in Eq.(18). 

In the same line, it is possible to write the expression of g{x) up to higher orders in r. For example, limiting the series 
expansion in Eq.(14) at third order in r, we found: 


gix) = 


rrieC 

kuT^. 


(ji - jo) + ^t(j 2 - 2ji -h jo) + - 3 j 2 + 3ji - jo) 


(28) 


Notice that while the expression derived for g(x) at first order approximation in r is independent of r, the expression 
for g{x) at higher order approximation in t depends directly on r. This is even more the case for the exact expression 
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of g{x) given in the following. In fact, using the exact form of the function P{s) given in Eq.(22), it is possible to write 
the exact form of the function g{x) as: 


9{x) = 


rrieC 

ksTe 


r-\-oc 


io(xe ^)P{s)ds — io{x) 


(29) 


The expression of g(x) approximated at first order in r as given by Eq.(p^ is the one to compare directly with the 
expression of g{x) obtained from the Kompaneets (1957) equation, since both are evaluated under the assumption 
of single scattering suffered by a CMB photon against the IC electrons. Fig.0 shows how the function g(x) tends 
to g{x) for lower and lower IC gas temperatures Tg. This confirms that the distorted spectrum obtained from the 
Kompaneets equation is the non-relativistic limit of the exact spectrum. Notice that the function g{x) given in Eq.(p^ 
is the spectral shape of the SZ effect obtained in the exact calculation while the function g{x) refers to the 1®* order 
approximated case of a single, thermal, non-relativistic population of electrons. 



Fig. 1. The function g{x) (solid line) is compared with the function g{x) for thermal electron populations with ksTe = 10 
(dot-dashed), 5 (dashes), 3 (long dashes) and 1 (dotted) keV, respectively. 


It is worth to notice that while in the non-relativistic case it is possible to separate the spectral dependence of the 
effect [which is contained in the function g{x)] from the dependence on the cluster parameters [which are contained in 
Compton parameter see eqs. (1-3)], this is no longer valid in the relativistic case in which the function Ji depends 
itself also on the cluster parameters. Specifically, for a thermal electron distribution, Ji depends non-linearly from the 
electron temperature Tg through the function Pi{s). This means that, even at first order in r, the spectral shape g{x) 
of the SZ effect depends on the cluster parameters, and mainly from the electron pressure Pth- 

To evaluate the errors done by using the non-relativistic expression g(x) instead of the relativistic, correct function 
g{x) given in Eq.(|^) we calculate the fractional error e = |[( 7 (a;) — g{x)]/g{x)\ for thermal populations with ksTe = 
10, 5, 3, 2 and 1 keV and for three representative frequencies (see Table 1). The fractional errors in Table 1 tend to 
decrease systematically at each frequency with decreasing temperature. Note, however, that the error found in the high 
frequency region, x = 15, is much higher than the errors found at lower frequencies and produces uncertainty levels of 
^ 50% for Te ^ 5 keV. This indicates that the high frequency region of the SZ effect is more affected by the relativistic 
corrections and by multiple scattering effects. In Table 2 we show the fractional error e done when considering the exact 
calculation of the thermal SZ effect and those at first, second and third order approximations in t for two values of 
the optical depth, r = 10“^ and t = 10“^, as reported in the Table caption. Even for the highest cluster temperatures 
here considered, ~ 20 keV, the difference between the exact and approximated calculations is ^ 0.25% at the 
minimum of the SZ effect {x ^ 2.3) and is ^ 2.23% in the high-frequency tail (a; ^ 15), the two frequency ranges 
where the largest deviations are expected and could be measurable, in principle. For high-temperature clusters, the 
third-order approximated calculations of the SZ effect ensures a precision ^ 2% at any interesting frequency. 
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ksTg = 10 keV 

ksTg = 5 keV 

ksTg = 3 keV 

fcsTe = 2keV 

ksTg = 1 keV 

X = 2.3 

7.20 

3.68 

2.26 

1.57 

0.93 

X = 6.5 

14.81 

7.28 

4.33 

2.87 

1.40 

X = 15 

64.84 

45.99 

32.67 

23.87 

13.15 


Table 1. The fractional error e = |[(7(a;) — g{x)]/g{x)\ computed for thermal electron populations with fcsTe =10, 5, 3, 2, 1 
keV is shown for different values of the frequency x. The quantity £ is given in units of 10“^. 



First order 

Second order 

Third order 

ksTg = 20 keV 




T = 10"2 




X = 2.3 

0.31 

0.22 

0.22 

X = 6.5 

0.12 

0.27 

0.27 

X = 15 

2.23 

1.71 

1.71 

fcflTe = 20 keV 




r = 10"® 




X = 2.3 

0.24 

0.23 

0.23 

X = 6.5 

0.21 

0.23 

0.23 

X = 15 

0.23 

0.18 

0.18 

fcflTe = 10 keV 




r = 10"2 




X = 2.3 

0.02 

0.02 

0.02 

X = 6.5 

0.7 

0.02 

0.02 

X = 15 

0.74 

0.21 

0.21 

fcflTe = 10 keV 




r = 10"® 




X = 2.3 

0.01 

0.01 

0.01 

X = 6.5 

0.01 

0.02 

0.01 

X = 15 

0.31 

0.25 

0.25 


Table 2. The fractional error e done considering the hrst, second and third order approximation in r compared with the exact 
calculation of the thermal SZ effect is reported for three interesting values of the frequency x. The assumed values of fcsTe and 
r are shown in the first column of the Table. The values of e are given in units of 10“^. 


Using the general, exact relativistic approach discussed in this section, we evaluate the frequency location of the 
zero of the thermal SZ effect, xq, for different cluster temperatures in the range 2 — 20 keV. The frequency location 
of Xq depends on the IC gas temperature (or more generally on the IC gas pressure) and in Fig. we compare the 
temperature dependence of xq evaluated at first order in r (which actually does not depend on r) with that evaluated 
in the exact approach for values r = 10“^ and t = 10“^. The location of the null of the SZ effect increases by 
~ 4.4% for ksTe up to ~ 20 keV. Assuming a quadratic fit Xq ~ a + bOg + where Qg = ksTg/vrigC^, we calculate 
the coefficients a,b and c which fit the temperature dependence of Xq. These are given in Table In particular the 



a 

h 

c 

Exact T = 0.001 

3.8271 

4.6059 

-4.8221 

Exact T = 0.01 

3.8262 

4.6704 

-5.4895 

Eirst order 

3.8294 

4.4304 

-1.7805 


Table 3. The coefficients of the quadratic fit of the behaviour Xo^th as a function of 9e 


ksTg/rriet? (see text for details). 


coefficients of the last row in Table ^ can be compared with those given by Itoh et al. (1998). At first order in ksTg, 
and for 1C gas temperatures up to 50 keV, these last authors found Xq ~ 3.830 + 4.4710e — 3.2680g. The coefficients 
given in the first two rows of Table ^ can be also compared with those found by Dolgov et al. (2000), who gave the 
expression of xq evaluated in an exact way for r, and for a temperature range 0 < ksTg < 50 keV and 0 < r < 0.05: 
xo = a(Tg) + Tl3{Tg) with a{Tg) = 3.830(1 + 1.1626le - 0.81446»2) and P{Tg) = 3.02106 - 8.67202. 

For T = 0.01, the expression for xq derived by Dolgov et al. (2000) writes as 


xo = 3.830 + 4.4810e - 3.2060^, 


(30) 
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Fig. 2. The location of xo for a thermal electron population is shown as a function of ksT,,. The exact calculation for r = 10 ^ 
(solid line) and for r = 10“® (dashed line) are shown together with the approximated calculation to first order in r (dotted 
line). 


while for r = 0.001 one gets 

xq = 3.830 + 4.4536»e - . (31) 

The remaining differences with our coefficients given in Table are due to the different Tg and r ranges used for the 
quadratic fit. 

3. The SZ effect generated from a non-thermal electron population 

Using the same general analytical approach derived in Sect. 2 above we derive here the exact spectral features of the 
SZ effect produced by a single non-thermal population of electrons. Such an exact derivation has been not done so far. 
For such a population we consider here two different phenomenological spectra: i) a single power-law energy spectrum, 
n^ei = like that which is able to fit the radio-halo spectra of many clusters, and ii) a double power-law 

spectrum which is able to fit both the radio halo spectrum and the EUV and hard X-ray excess spectra observed in 
some nearby clusters (see Petrosian 2001). Our formalism is, nonetheless, so general that it can be applied to any 
electron distribution so far considered to fit both the observed radio-halo spectra and those of the EUV/Hard X-ray 
excesses (see, e.g., Sarazin 1999, Blasi & Colafrancesco 1999, Colafrancesco & Mele 2001, Petrosian 2001). 

The single power-law electron population is described by the momentum spectrum 

fe,rei{p;Pi,P 2 ,a) = A{pi,p 2 ,a)p~°‘ ; P1<P<P2 (32) 

where the normalization term A{pi,p 2 ,a) is given by: 

Mpi ,P2,a) = (33) 

Pi -P2 

with a ~ 2.5. This is the simplest electron distribution which is consistent with the spectra of the radio halos observed 
in many galaxy clusters (see, e.g., Feretti 2001 for a recent review). In the calculation of the non-thermal SZ effect 
we consider the minimum momentum, pi, of the electron distribution as a free parameter since it is not constrained 
by the available observations, while the specihc value of p 2 3> 1 is irrelevant for power-law indices a > 2 which are 
required by the radio halo spectra observed in galaxy clusters. 

The distribution in Eq.(^) can be however affected by a crucial problem: if the electron spectrum is extended at 
energies below 20 MeV (i.e., at p ^ 40) with the same spectral slope a ~ 2.5 of the high energy tail, the heating rate of 
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the IC gas produced by Coulomb collisions of the relativistic electrons becomes larger than the bremsstrahlung cooling 
rate of the IC gas (see, e.g., Petrosian 2001). This would produce unreasonably and unacceptably large heating of the 
IC gas. 

For this reason, we consider also a double power-law electron spectrum which has a flatter slope below a critical value 
Per- 


fe^rel(,P: Pl^ P2 : Per ; ttl, O 2 ) ^ {pi , P2 7 Per: , 0^2 ) 


P Pi <P < Per 

p--ai+a2 p-a2 p^^ < p <_ p^ 


(34) 


with a normalization factor given by: 


K{pi,p2,per7 ai, 02 ) = 


1—ai 

Pi 


-p. 


1—ai 


-OL\-\-cx.2 


1 — 0:2 


1 - 

P 2 


ai — 1 


02 — 1 


-1 


(35) 


If Per ^ 400, the electron distribution of Eq.(|3j) with ai ~ 0.5 at p ^ Per can be extended down to very low energies 
without violating any constraint set by the IC gas heating. Hence, we assume here the following parameter values: 
ai — 0.5, Q !2 = 2.5, Per = 400 e p 2 —> oo. We again consider pi as a free parameter. 

A crucial quantity in the calculation of the non-thermal SZ effect is given by the number density of relativistic 
electrons, rie^rei- The quantity ne,rei in galaxy clusters can be estimated from the radio halo spectrum intensity, but 
this estimate depends on the assumed value of the IC magnetic field and on the model for the evolution of the 
radio-halo spectrum (see, e.g. Sarazin 1999 for time-dependent models and Blasi & Colafrancesco 1999, Colafrancesco 
& Mele 2001 for stationary models). For the sake of illustration, the radio halo flux J^, for a power-law spectrum 
is, in fact, given by Jj/ oc ne,reiB°‘^~^v~°‘^, where ar is the radio halo spectral slope. Because of the large intrinsic 
uncertainties existing in the density of the relativistic electrons which produce the cluster non-thermal phenomena, 
a value Ue^rei = 10“® cm“^ for pi = 100 has been assumed here. Our results on the amplitude of the non-thermal 
SZ effect can be easily rescaled to different values of ne,rei- in fact, decreasing (increasing) ne,rei will produce smaller 
(larger) amplitudes of the SZ effect, as can be seen from eqs.(22,24,28,29). 

The density Ug^rei increases for decreasing values ofpi. In fact, multiplying the electron distribution in Eq.(^) by the 
quantity Ue^reiipi) one obtains 


Ne{p-,Pl) = ne,rel{Pl)A{Pl)P “ 

where the function A{pi) is given by Eq.(^. Thus, the electron density scales as 


n'e,re/(Pl) — '^e,re[{pi) 


Ajpi) 

A{pi) 


(36) 


(37) 


where we normalized the electron density at a fixed value of pi = 100. 

The value of pi sets the value of the electron density as well as the value of the other relevant quantities which 
depend on it, namely the optical depth Trei and the pressure Prei of the non-thermal population. In particular, the 
pressure Prei for the case of an electron distribution in Eq.(|2|), is given by 


Prel 


dpfe{p)^Pv{p)meC 


UemeC^ia — 1 ) 


6[p 


l-alPi 
P2 


B_ 


a — 2 3 — a 


Pi 

P2 


(see, e.g., Ensslin & Kaiser 2000), where is the incomplete Beta function 


(38) 


Bx{a,b) = f P — ^dt 

Jo 


(39) 


(see, e.g., Abramowitz & Stegun 1965). The optical depth of the non-thermal electron population is given by 


Trel{Pl) =2-10 


6 Mp'i) 

ne,rel{Pl = 100) 

e 

A{pi) 

10 “® cm~^ 

1 Mpe 


(40) 


For an electron population as in Eq.(^) we obtain analogous results. The electron density with the same normalization 
as before is given by 


^e,re/(Pl) — ^e,rez(Pl) 


KjPi) 

K(.Pi)' 


(41) 
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-50 5 10 15 20 25 


s 


Fig. 3. The distribution Pi(s) for a non-thermal population as in Eq.(p^ with a = 2.5, p 2 ^ oo and pi = 0.5, 1, 10, 100, 1000, 
as indicated. 


The optical depth is again given by Eq.(40), where we substitute the ratio K(pi)/K(pi) to A{'pi) / A(pi) . The electron 
pressure of this last population has instead a more complicated expression given by: 


Prelipi) — 


K (j)i)ne{pi)meC^ 


6 


B 


B 


a2 — 2 3 — 02 


ai — 2 3 — 01 




-I Pc 


(42) 


To evaluate the SZ effect induced by a single non-thermal electron population we use the relativistically correct 
formalism described in Sect. 2 which also takes into account the effects of multiple scattering. Again, we are working 
here in the Thomson limit. The general expression for the distorted spectrum I{x) is given again by Eq.( 8 ). The only 
change in our general formulation consists in inserting the non-thermal distribution fe,non-th in the function T’i(s). 
The functions Pn{s), which determine the probabilities to have a frequency shift s due to n scattering are still given 
by the convolution of n distributions Pi{s) (see Eq. 6 ). Fig.|| shows the function Pi{s) for different values of pi in the 
range 0.5 — 1000. For low values of pi, the function Pi{s) is centered around s = 0 and does not produce, hence, 
large frequency changes. For higher values of pi (i.e., for increasing values of the minimum energy of the electrons) 
the function Pi{s) becomes wider and it is centered at higher and higher values of s. For pi = 1000 the redistribution 
function Pi{s) is centered around s = 15, which corresponds to a frequency change v'jv = ^ Z - 10®. Thus, the SZ 

effect caused by electronic populations with very high values of pi moves the photon frequencies from the CMB region 
up to very different regions causing a depletion of scattered CMB photons at sub-mm wavelengths where the SZ effect 
is usually studied. We report in Table ^ the values of the pressure and of the density of a single power-law population 
as a function of pi. 

To emphasize the changes in the photon redistribution function produced by a non-thermal electron distribution, 
we show in Fig.^ the functions Pi{s),P2{s) and P-i{s) for a single and double power-law electron distributions. It is 
evident that these functions are wider and more a-symmetrically skewed towards large and positive values of s with 
respect to the same distributions evaluated in the case of a single thermal electron population. 

Once the exact redistribution function P{s) is known, it is possible to evaluate the exact spectral distortion AI{x) 
through Eq.(||). 

It is also possible to evaluate the distorted spectrum using Eq.(16) which gives the series expansion up to the n-th 
order in r and to compare the exact and approximated shapes for I{x). We show in Figj^the first three terms &i, 62 , ^3 
of the series expansion in Eq.(16) as a function of the a-dimensional frequency x. The approximated expression for 
/ (a:) up to third order in t is an excellent fit of its exact behaviour at frequencies around the minimum of the effect 
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Pi 

Prei (keV cm""*) 

ne.rei (cm“"') 

0.5 

0.59 

2.83 • 10“^ 

1 

0.47 

1 •10"® 

10 

0.16 

3.16 • 10~® 

100 

5.11 • 10"^ 

1•10"® 

1000 

1.61 • 10"^ 

3.16 • 10“® 


Table 4. Values of the pressure and of the density of a single power-law population as in Eq.(|3^ with a = 2.5, p 2 —> oo and 
ne,rei{pi = 100) = 1 • 10~® cui”® for different values of pi. 




S S 


Fig. 4. The distributions Pi{s) (solid line), P 2 (£l (dashed) and P 3 {s) (dotted) evaluated for non-thermal electron populations 
with a single power-la w sp ectrum given in Eq.(p^ with a = 2, pi = 1 and p 2 = 100 (left panel) and with a double power-law 
spectrum given in Eq.(M) (right panel). Note that these distributions are wider and more skewed towards larger values of s 
than the corresponding distributions evaluated in the case of a thermal electron population. 


(x ~ 2.3) and it shows a difference ~ 3% at a; « 6.5 and a difference « 0.1% at high frequencies, a: ~ 15 (see TableH). 
The precision of the approximation depends also on the value of r and we verified that the precision of the second and 
third order approximations does not increase for values r ^ 10“^ at any frequency between a; ~ 2 and x ~ 15. 



First order 

Second order 

Third order 

r = 0.01 




X = 2.3 

0.82 

0.44 

0.44 

X = 6.5 

1.42 

2.90 

2.90 

X = 15 

0.83 

0.12 

0.12 

r = 0.001 




X = 2.3 

0.49 

0.45 

0.45 

X = 6.5 

2.73 

2.88 

2.88 

X = 15 

0.19 

0.12 

0.12 


Table 5. Fractional difference (in units of 10 ^) between the exact spectral distortion and that evaluated at first, second and 
third order in r. We report here the cases for values of the optical depth r = 0.01 and r = 0.001. 
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Fig. 5. The functions bi{x) (solid), 62 ( 2 :) (dashed) and 63 ( 2 ;) (dotted) evaluated for a non-thermal electron population with a 
single power-law spectrum as in Eq.(32). 


3.1. The function g{x) for a single non-thermal population 


Also for a non-thermal electron population it is possible to write the spectral distortion in the general form of Eq.(22) 
as: 




-th 


ix)=2 


{keTof 

[hcf 


Vnon—th ; 


where the Comptonization parameter is given by the general expression 


Vnon — th — 


(JT 

rrieC^ 




(43) 


(44) 


in terms of the electronic pressure Prei which, for a single and double power-law populations, is given by Eq.(^) and 
by Eq.(42), respectively. 

At first order in t we can write the spectral function g{x) of the non-thermal SZ effect as 


- jo] _ . 

y ^ J PM - 

where we defined the quantity 

jpde 


(45) 


(fcsTe) = ^ J Pd£ 


J riedi 


= J dpfe{p)-pv{jp)meC , (46) 

which is the analogous of the average temperature for a thermal population. In fact, comparing Eq.(^^ with Eq.(^^, 
{ksTg) = ksTe obtains. For a non-thermal population with a single power-law distribution one gets 


(fcsTe) = 


meC^{a — 1 ) 
6[p 


l-alPi 
\P2 


B 


-2 3-a 


(47) 


while in the case of a double power-law distribution we find instead 


(fcsTe) = 


K{p 


i)me(r 


6 


B 


Ofi — 2 3 — Qi 


+ 


-J Pc 
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B 



— 2 3 — Q!2 
2 ’ 2 


-OL\-\-CX.2 


-I P2 


(48) 


It is possible to write the Comptonization parameter ynon-th as a function of the quantity {ksTe) and of the optical 
depth r: 


(Tt 

Unon—th — o 

rrieC^ 


Preldi = (Jt 


{ksTe) 

nieC^ 



{ksTe) 

2 ^ 
nieC^ 


(49) 


(note that here r = r^ez holds). 

We can also write the n-th order approximations of the distorted spectrum using Eq.(16). For example, at third order 
in T, the following explicit expression founds: 


777 / 1 1 

= TTtVt (-^1 “ - 2ji + jo) + - 3j2 + 3ji - jo) 

{ksTe) [2 6 


(50) 


Finally, the exact form of the function g{x) for a non-thermal electron population writes as: 


Kx) 


rUeC^ (1 
(ksTe) {t 


^+oo 


io{xe '^)P{s)ds — io{x) 


(51) 


We show in Fig.^a comparison between the function g(x) obtained under the Kompaneets approximation [see Fq.(2)] 
and the function g{x), approximated to first order in r, and evaluated for a single power-law population. A major 
difference between the two functions is the different position of the zero of the SZ effect which is moved to higher 
frequencies in the case of the non-thermal population with respect to the case of the thermal one. In Fig.^we compare 
the function g(x) for a non-thermal population evaluated at first order in r (which is actually independent from the 
value of r) with the exact function for values r = 1, r = 0.1 and r = 0.01. We notice that the first order approximation 
is a good approximation of the exact case for values r < 10“^. Fig.|| shows the function g{x) (evaluated at first order 



X 


Fig. 6. The spectral function g{x) for a thermal population in the Kompaneets limit (solid line) is compared with the function 
g{x) approximated at hrst order in r for a non-thermal pop ulation with a single power-law spectrum with parameters pi = 4, 
P 2 = 10^° and a = 2.5 (dashed line), evaluated from Eq.(|45[). The dotted line represents the zero reference value. 


in r) for a single power-law population and considering different values of pi. For increasing values of pi the maximum 
of the function g{x) moves towards higher and higher frequencies. As a consequence, also the zero of the non-thermal 
SZ effect moves to higher values of x as shown also in Fig.^j. The location of xo changes rapidly with increasing values 
of Pi and reaches a value xo « 22 for pi = 2 - 10^. At such values of pi the pressure Prei takes values ~ 0.02 keV cm“^. 
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X 

Fig. 7. The function g{x) for a non-thermal population with single power-law spectrum with parameters pi = 4, p 2 = and 
a = 2.5. We show the first order approximation in r (solid line) and the exact case for r = 1 (dashed line), r = 0.1 (dotted 
line) and r = 0.01 (dot-dashed line). Note that the first order approximation is a reasonable approximation of the exact case 
only for small values of r. 




Fig. 8. In the left panel we show the function g{x) (see Eq.^^ for a non-thermal population as in Eq.(^^ with pi = 0.5 (solid 
line), 1 (dashed) and 10 (dotted). The right panel shows the same function for p\ = 10 (solid line), 100 (dashed) and 1000 
(dotted). 


as shown in Fig.9. An important point to notice is that the electron density of the non-thermal population decreases 
for increasing values of pi as was already shown in Table For this reason the pressure Prei and the resulting spectral 
distortion decreases in amplitude for increasing values of pi. This result is shown more specifically in Figj^ where we 
plot the quantity /S.i{x) = A/(x)[(hc)^/2(fc_BTb)^] for values of pi in the range 0.5 — 10^. 
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Pi P 


Fig. 9. The behaviour of the zero of the SZ effect for a non-thermal population given in Eq.(|3^ as a function of pi (left panel) 
and of the pressure Prei (right panel) expressed in keV cm“®. Note that the pressure Prei decreases for increasing values of pi, 
as discussed in Sect.3 



Fig. 10. The spectral distortion Ai{x) in units of 2{kBTo)^/{hc)^ for a non-thermal population as in Eq.(p2|) evaluated for 
Pi = 0.5 (solid line), 1 (dashes), 10 (dotted), 100 (long dashes) and 1000 (dot-dashes). 


For an electron population with a double power-law spectrum as in Eq.(|^) the values of the density and pressure 
are also reported in Table || as a function of pi (note that the case pi = 1000 is the same of the single power-law 
spectrum). The values of Prei and ne,rei for low values of pi are now smaller than those of the single power-law case 
(see Tabled), thus avoiding the problem of an excessive heating of the cluster IC gas for low values of pi. 

In Fig.^ we show the function Pi{s) evaluated for an electron population with double power-law spectrum (see 
Eq.34) and its dependence on the lower momentum cutoff pi. The behaviour of Pi{s) in this case is quite different with 
respect to the single power-law case (see Fig.||) and depends mainly on the flatness of the electron spectrum below the 
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Pi 

Prei (keV cm~'=‘) 

We.rei (cm~'=‘) 

0.5 

1 

10 

100 

1.26 • 10"^ 

1.26 • 10"® 

1.23 • 10"® 
9.88- 10-'‘ 

1.56 • lO"*’ 
1.54 • 10"® 
1.41 • 10"® 

1 •10“® 


Table 6. Values of the pressure and of the density of a double power-law population as in Eq.(^) with ai = 0.5, Qf 2 = 2.5, 
Per = 400, P 2 —> C )0 and Ue^reiipi = 100) = 1 • 10“® cni“® for different values of pi. 



-50 5 10 15 20 25 


s 


Fig. 11. The function Pi{s) [see Eq.(^)] is shown for non-thermal populations given in Eq.(^^ with pi = 0.5 (solid), 1 (dashes), 
10 (dotted), 100 (dot-dashes). 


momentum cutoff Per- We notice that i) even for low values of ~ 0.5 small frequency changes are unlikely, as shown 
by the strongly a-symmetric shape of Pi{s); ii) there are small differences in the function Pi{s) for different values of 
Pi ~ Per because the electron density in this range is lower than in the single power-law case and its variation is mild. 
We also show in Fig.|^ the function g{x) evaluated at first order in r for a double power-law electron population as 
in Eq.(^^. In this case, there is no well defined maximum of the SZ effect since low frequency changes are unlikely. 
Moreover, the function g{x) is less sensitive to different values of the lower momentum cutoff pi for pi < per because 
of the flatness of the electron spectrum in this region. This fact is also confirmed by the variation of the zero location 
of the non-thermal SZ effect (see Fig.|Tl). The location of xq changes now less rapidly than in the case of a single 
power-law distribution (see Fig.11) due to the different shape of g{x), and hence of P(s), caused by the lower density 
of the low-energy electrons with respect to the case of the single power-law spectrum. The zero of the non-thermal 
SZ effect produced by a single non-thermal population is also shifted to frequencies Xg ~ 10 (much higher than the 
value Xo,th = 3.83) due to the large frequency shifts experienced by the CMB photons scattering the high energy 
non-thermal electrons. 

We finally show in Fig.|^ the spectral distortion Ai{x) for different values of pi. Note, finally, that the amplitude of 
AI{x) produced by an electron population with a double power-law spectrum as in Eq.(^) is much lower than that 
produced by a population with a single power-law spectrum as in Eq. (^^. 

The case of a system dominated by non-thermal electrons is not applicable straightforwardly to galaxy clusters but 
rather to radio galaxies whose jets inject in the ICM large quantities of relativistic non-thermal electrons. We want to 
stress that our general approach is able to describe the Compton scattering distortions even in the cases of the radio 
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Fig. 12. The function g(x) [see Eq.(|45|)] for non-thermal populations as in Eq.(^) with pi = 0.5 (solid line), 1 (dashes), 10 
(dotted) and 100 (dot-dashes). 



Fig. 13. The behaviour of the zero of the SZ effect for a non-thermal population with spectrum given by Eq.(^) as a function 
of Pi (left panel) and of the pressure Prei (right panel) in units of keV cm“^. As in Fig.9, the pressure Prei decreases for 
increasing values of pi but with a different rate depending on the different spectrum of the electron distribution. 


lobes and of the jets of radio galaxies or AGNs with large optical depths and ultra-relativistic electron energies. We 
will address more specifically this issue in a forthcoming paper. 

4. The total SZ effect produced by a combination of thermal and non-thermal populations 

For galaxy clusters which contain two (or more) different electronic populations, like radio-halo clusters which contain 
a thermal population emitting X-rays and a non-thermal population producing the radio halo emission, one has to 
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Fig. 14. The spectral distortion in units of 2(fcsTo)®/(/ic)^ for a non-thermal population given by Eq.(^) computed for p\ = 
0.5 (solid), 1 (dashes), 10 (dotted) and 100 (dot-dashes). 


evaluate the spectral distortion produced simultaneously by both populations on the CMB radiation. Such a derivation 
has been not done so far. We assume here that the two populations are independent and that no change on the 
thermal population is induced by the non-thermal electrons. This condition is reasonable for non-thermal electrons 
with energies 2 150 MeV which do not appreciably loose energy through Coulomb collisions in the ICM (see, e.g., 
Blasi & Colafrancesco 1999; Colafrancesco & Mele 2001) and hence do not produce a substantial heating of the IC 
gas. Electrons with Eg 2 150 MeV loose energy mainly through ICS and synchrotron losses. The radio halo emission 
in galaxy clusters is, in fact, produced by electrons with energy E,, k 16.4 GeV {BjjGH (see, e.g., 
Colafrancesco & Mele 2001) which yield Eg ^ 1.6 — 52 GeV for a typical IC magnetic field at the level oi B = 1/iG 
at the typical frequencies, Vr ^ 10“^ — 10 GHz, at which radio halo spectra are observed. For such energies the 
non-thermal electrons do not strongly affect the thermal IC gas. However, we discuss in the following the impact of 
different non-thermal population spectra on both the total SZ effect and on the cluster IC gas. 

Under the hypothesis of independence between the thermal and non-thermal electron populations, the probability 
that a CMB photon is scattered by an electron of population A (say thermal) is not affected by the fact that the 
same photon has been scattered by an electron of population B (say non-thermal). Thus, at first order in r we expect 
that the total SZ effect is given by the sum of the first order SZ effects due to each single population. At higher order 
approximation in r, however, one has to consider the effect of repeated scattering (see Sect.2). It is intuitive to think 
that in this last case the total SZ effect is not simply given by the sum of the separate SZ effects at higher orders in r 
and the terms describing the cross-scattering between CMB photons and the electrons of populations A and B have 
to be taken into account. In the following, we will derive in details the general expression for the total SZ distortion 
caused by a combination of two electronic populations. 

Let us consider two electron populations with momentum distributions Ja (p) and Jb (p) with optical depths ta and 
tb, respectively. We also assume that the two electron distributions are both normalized to 1. The total distribution 
of the electron momenta is 

fe{p) = CAfAip) + CBfsip) , (52) 

with appropriate normalization coefficients ca and Cb- Since the total distribution fe{p) is also normalized to 1, the 
coefficients ca and cb are related by the condition -I- cb = I- The coefficients ca and cb yield the weight of each 
population with respect to the other and thus their ratio is given by 

CA _ Ta 
Cb Tb 


(53) 
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Such a condition, together with the normalization condition, yields: 
ta tb 

CA = -;- CB = -;-• 

ta + tb ta + Tb 

The total distribution fe{p) has an optical depth given by 


(54) 


T = ta +Tb. 


(55) 


In fact, the probability that a CMB photon can suffer n scatterings by the electrons of each distribution is given by the 
sum of all the possible combinations of the probability to suffer ua scatterings from the electrons of the distribution 
Ja and the probability to suffer ub scattering from the electrons of the distribution /s, with ua + ns = n^. Since the 
two distributions are independent, this probability is given by the product of the two separate probabilities: 


Pn = 


nA+nB=n 


E 

L+raS 

E 

L+' 

n 

E 


VriA ’ Vni 


= E 


nA+nB=n 


ua'- 


Ub'- 


„-{TA+TB)~nA UB 


nA+nB=n 


n^!ns! 


nA=0 


^-(ta+tb) 


„-(TA + TB)~nA n-UA 

UA^in - HaY- 


E 


ni 


^ nl nAl{n - n^)! ^ 

riA—O ^ ' 


e + tbY 


n\ 


n\ 


(56) 


(we assume here a Poisson probability distribution). This means that the electrons of the distribution /b simply add 
to the electrons of the distribution /a- In fact, from the definition of optical depth, one gets: 


T = <7t I Uedi 


= TA + Tb = (Tt jiriA + nB)d£. 


(57) 


The exact expression for the spectral distortion produced by two electron populations can be evaluated from eqs. 
-(21), using the appropriate expression for Pi{s). The general expression for the distribution Pi{s) is given by: 


POO 

Pi{s) = / dpfe{p)Ps{s;p). 
Jo 


(58) 


Inserting Eq.(|52|) in Eq.(p8|) one obtains the function Pi{s) for two populations as a function of the distributions Pi{s) 
of each single population: 


pOO pOO 

Pi{s) = CA dpfA{p)Ps{s;p)+CB / dpfB{p)Ps{s;p) 

Jo Jo 

= c^-Pia(s)-I-cbPi_b(s). 

The Fourier transform of Eq. (^^ is 

Pi{k) = CAplAik) + CBplB{k) 

and from this expression one obtains the Fourier transform of the total redistribution function: 

P{k) = = 

= ^-Pa+tb) [l- ^YY-tb pA(k)- ^YY tb Pig (^)] = 

_ g-'rA[i-PlA{k)]-TB[l-PiB{k)] _ 

= pA{k)pB{k). 


(59) 

(60) 


(61) 
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The exact total redistribution function Ptot{s) is given by the convolution of the redistribution functions of the separate 
electron populations: 

Ptot{s) = Pa{s) 0 Pb{s). (62) 

Thus, the exact spectral distortion produced by two electron populations on the CMB radiation is given by: 


/ + 00 

Io{xe~'‘)Ptot{s)ds, 

-CXD 


(63) 


in terms of the exact total redistribution function Ptot{s) given in Eq.(62). We also give in the Appendix the analytic 
expressions for the approximations at first and second order in r of the total spectral distortion Itot{x). 


4.1. The function g{x) for a combination of two electron populations 

The total SZ effect produced by the combination of two electron populations can be written again in terms of a 
Comptonization parameter ytot and of a spectral function g{x), in the form of Eq.(22). For two electron populations, 
we can write the average temperature in Eq.(^^ as: 


/.CX) ^ 

(fcsTe) = J dpfe{p)-pv{p)meC = 


1 


= J dp[cAfAip) +CBfB{p)]-^pv{p)meC = 

= CA{kBTe)A+CB{kBTe)B- 

The Comptonization parameter ytot can be evaluated using Eq.(^) as: 
1 


(64) 


ytot — 


TTleC^ 

1 


:[cA{kBTe)A + CB{kBTe)B]{'rA + Tb) = 
[TA{kBTe)A + TB{kBTe)B\ = 


rrioC^ 

= yA + VB- 

At Hrst order in r, the function g(x) for two populations can be expressed as 
^i{x) TA[jiA{x) - jo(a;)] + TB[jiB{x) - jo{x)] 


(65) 


gix) = 


ytot 

VAgAjx) + yBgBjx) 

VA + VB 


VA + VB 


( 66 ) 


in terms of the functions g{x) for each single electron population. We notice that in this case the function g{x) depends 
from the Comptonization parameters, and hence from the optical depths, of the two electron populations even at the 
first order approximation in r. Moreover, at higher orders in r, the function g{x) cannot be expressed simply in terms 
of the spectral functions of the single populations because in the expression for Ai{x) [see Eq. ( A.10| )], cross-correlation 
terms appear. 


4.2. The SZ effect produced by a combination of a thermal plus a non-thermal populations 

Here we apply the previous derivation of the total SZ effect to a galaxy cluster which contains both a termal and a 
non-thermal electron populations. In the following computation of the total SZ effect for a Coma-like cluster we choose 
a thermal population with the following parameters: fcsTe = 8.5 keV, neAh ~ 3 • 10 


^ cm ^ and a cluster radius of 


1 h^Q Mpc. With these parameters the optical depth of the thermal population is 

Tth-b - 10” 


3 

'^e,th 


i 


.3 • 10“^ cm~^. 


AKo Mpc. 


and its pressure is: 

Pth ^ 2.55- 10-2 fceCcTO- 


3 

’^e^th 


kBTe 


-3 • 10“^ cm~^. 


8.5 keV 


(67) 


( 68 ) 


As for the spectrum of the non-thermal population we consider both the phenomenological cases of a single and double 
power-law populations [see Eq.(p^ and Eq.(p^ respectively). The parameters we consider are specifically: a = 2.5 
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Fig. 15. The function Pi{s) [see Eq. (^)] for a combination of thermal population and non-thermal population as in Eq.(^^ 
with Pi = 0.5 (solid), 1 (dashes) and 10 (dotted). 


and p 2 ^ oo for the single power-law population and ai = 0.5, 02 = 2.5, with pcr = 400 and p 2 ^ 00 for the double 
power-law population. We consider pi as a free parameter in both cases. The relativistic electron density has been 
normalized at ne,rei{pi) = 10“® cm“^ for pi = 100 in both cases. The final SZ effect can be re-scaled simply to different 
values of Ue^rei as discussed in Sect. 3. 

The amplitude of the non-thermal the contribution to the total SZ effect increases with decreasing value of pi. For 
the case of a single non-thermal power-law population it becomes appreciable for pi 2 10 while for pi 2 100 the total 
SZ effect is indistinguishable from the pure thermal effect. In fact, for pi 2 10 the function Pi{s) is nearly coincident 
with that of a single thermal population as shown in Fig.|l5|, while for pi < 10 the function Pi(s) is more a-symmetric 
towards large and positive values of s. The peaks of the non-thermal distributions at large s are present in the total 
distribution Pi{s) but their amplitude decreases with Trei and becomes lower and lower with increasing values of pi. 

The amplitude and the shape of the total SZ effect in a cluster which contains non-thermal phenomena depends 
on the spectral and spatial distribution of the non-thermal population, given that the cluster IC gas properties are 
well known from X-ray observations. We can take advantage of this fact to use the non-thermal SZ effect as a tool to 
constrain the physical properties of the non-thermal population. We show in Fig.16 the total spectral distortion due 
to the SZ effect in a cluster with non-thermal population as a function of pi. Low values ofpi 2 1 for the non-thermal 
population induce large distortions with respect to the pure thermal SZ effect. However, for low values of pi the 
pressure of the non-thermal population with single power-law distribution exceeds substantially the thermal pressure 
Pth (see Table |^) ; as a consequence, the low energy electrons produce an exceedingly large pressure unbalance and 
also a large heating of the IC gas, which are not acceptable. 


pi 

Prel /Pth 

'^e,th{'^e,re.l 

0.5 

23.3 

1.06 

1 

18.6 

3.00 

10 

6.33 

94.87 

100 

2.00 

3000 

1000 

0.63 

94868 


Table 7. Values of pressure and density ratios for a combination of thermal population and a non-thermal population with 
single power-law spectrum for different values of pi (see text for details). 


The frequency location of the zero of the total SZ effect, xq, is also a powerful diagnostic of the presence and of 
the nature of the non-thermal population. The frequency location of Xq increases for decreasing values of pi as the 
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Fig. 16. The spectral distortion, in units of 2(fc_BTo)®/(/ic)^, of a single th erma l population of a Coma-like cluster (solid) is 
compared with the spectral distortion evaluated at hrst order in r (see Eq. (A.6)), produced by a combination of thermal and 
non-thermal population given in Eq.(^^ with pi = 0.5 (dashes) and 1 (dotted) (left panel) and for pi = 10 (dashes), 100 
(dotted) and 1000 (dot-dashes) (right panel). 




Fig. 17. The behaviour of the zero of the total SZ effect evaluated for a combination of thermal and non-thermal population 
given by Eq.(^2|) is shown as a function of minimum momentum, pi, of the non-thermal population (left) and of the pressure 
ratio Prei/Pth (right). 


impact of the non-thermal population becomes larger (see Fig.|l^). Consistently, Xo increases with increasing values of 
the pressure ratio Prei/Pthj as also shown in Fig.|^. 

In the case of a double power-law non-thermal population the ratio Prei/Pth always remains of order 10“^ even for 
vary low values pi ^ 0.5 (see Table ^), and in this case there is a negligible dynamical and thermal influence of the non- 
thermal population on the thermal one. In this case, also the total distribution f’i(s) and the total spectral distortion 
are little affected by the value of pi. We show in FigjT| the spectral distortion of a single thermal population compared 
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Pi 

Prel/Pth 

'^e,th / '^e,rel 

0.5 

4.94 • lO”"" 

1926 

1 

4.94 • 10“^ 

1948 

10 

4.84 • 10“^ 

2127 

100 

3.87-10-2 

3000 


Table 8. Values of pressure and density ratio for a thermal population and a non-thermal population with double power-law 
spectrum for different values of pi (see text for details). 


with the spectral distortion produced by a combination of thermal and non-thermal population with double power-law 
spectrum. The contribution of the non-thermal population to the total SZ effect depends on the normalization of the 
density ne,rei (which is highly uncertain). For increasing values of the density of the non-thermal population the non- 
thermal SZ effect becomes more relevant and produces substantial changes to the SZ distortion produced by a single 
thermal population at frequencies a; 10 (see Fig.|l^. Much milder changes are present in the high frequency range 
a; ^ 15 due to the low amplitude of the non-thermal signal in this region (see, e.g., Fig. 14). We also show in FigJ^ 


ne,rel{pi = 100) (cm-"‘) 

Prel/Pth 

! '^e,rel 

Xo 

1 •10-“ 

4.94 • 10-2 

1926 

3.9085 

1 •10-® 

0.49 

192 

3.9965 

3-10-® 

1.48 

64 

4.1765 


Table 9. Values of pressure aud density ratio and the corresponding position of the zero for a thermal population and a 
non-thermal double power-law population for pi = 0.5 as a function of ns^reiipi = 100) (see text for details). 


the behaviour of xq as a function of pi and of the pressure ratio for a density normalization ne,rei{pi = 100) = 3 • 10“® 
cm“^. High values of the ratio Prei/Pth produce a substantial displacement of xq up to values xq ^ 4 which correspond 
to frequencies ^ 227 GHz. In Fig.|T| we also show the behaviour of xq a a function of ne,rei{pi = 100) for pi = 0.5. 
We also report in Table § the values of the density and the corresponding values of Prei/Pth and of ne,thlne,rei as 
well as the frequency position of the zero of the total SZ effect. It is interesting to note in this context that the 
precise detection of the frequency location of xq provides a direct measure of the ratio Prei / Pth and, in turn, of the 
energy density and of the properties of the spectrum of the non-thermal population once the pressure of the thermal 
population is known from X-ray observations. In fact, all these physical information contributes to set the specific 
amplitude and the spectral shape of the total SZ effect (see Figs.17 and 20 and Figs. 16, 18 and 20 for the different SZ 
effect features produced by two different non-thermal electron populations). This is a relevant and unique feature of the 
non-thermal SZ effect in clusters since it yields a measure of the total pressure in relativistic non-thermal particles in 
the cluster atmosphere, an information which is not easily accessible from the study of other non-thermal phenomena 
like radio halos and/or EUV or hard X-ray excesses. Thus, the detailed observations of the frequency shift of xq in 
clusters provides unambiguously relevant constraints to the relativistic particle content in the ICM and to their energy 
distribution. Such a measurement is also crucial to determine the true amount of kinematic SZ effect due to the cluster 
peculiar velocity since it is usually estimated from the residual SZ signal at the location of the zero of the thermal 
(relativistic) SZ effect. Due to the steepness of the SZ spectral shape in the region of the null of the overall SZ effect, 
the possible additional non-thermal SZ signal must be determined precisely in order to derive reliable limits to the 
kinematic SZ effect. We will present a more detailed analysis of the kinematic SZ effect in a forthcoming paper. 

5. The SZ effect produced by a combination of different thermal populations. 

The general derivation of the total SZ effect produced by different electron populations described in Sect. 4 above 
can be generally applied also to the case of a combination of different thermal populations. So, in this Section we 
consider the SZ effect produced by two thermal populations with different temperatures and densities residing in 
the same cluster. The presence of an additional thermal IC population of quite low temperature, T ^ 1 keV, has 
been considered as a possibility to explain the origin of the EUV/soft X-ray excess observed in some nearby clusters 
(Coma and Virgo for instance, see e.g. Lieu et al. 1999) and it is also expected in the outer regions of the clusters 
due to infall of warm gas from large scale filaments. Here we want to explore the modifications that such additional 
warm component can introduce into the overall SZ effect in galaxy clusters. More interestingly, we want to show that 
is possible to obtain constraints on any additional thermal component in galaxy clusters using SZ observations. We 
consider for this additional thermal population temperatures in the range kBTe ^2 = 0.5—1 keV and densities ne ,2 in the 
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Fig. 18. The spectral distortion, in units of 2{kBTo)^/{hc)^, of a single thermal population (solid) is compared with the spectral 
distortion evaluated at first order in t produced by a combination of thermal and non-thermal population given in Eq.(p4|) with 
Ue.reilPi = 100) ” 1 ' 10“® and Pi = 0.5 (dashes), 1 (dotted), 10 (dot-dashes) and 100 (long dashes). We show the enlargement 
in the regions of the minimum (top), of the zero (mid) and of the maximum (bottom) of the total SZ effect. 
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Fig. 19. The spectral distortion, in units of 2{kBTo)^/{hc)^, of a single thermal population (solid) is compared with the spectral 
distortion evaluated at first order in r produced by a combination of thermal and non-thermal population given in Eq.(p^ with 
Pi = 0.5 and ne{pi = 100) = 10“® (dashes) and 3 • 10“® cm“® (dotted). 


range 10“^ — 10“^ cm“^. We also assume here, for simplicity, that the two thermal populations have the same spatial 
distribution. Moreover, we report in Tables and ^ the pressure ratio, the density ratio and the frequency location 
of the zero of the total SZ effect for two values of the temperature of the warm population: kBTe ^2 = 0.5 and 1 keV, 
respectively. The presence of additional cooler electrons produces a tightening of the photon re-distribution function 


ne ,2 (cm "*) 

P 2 /P 1 

Ue,l/ne,2 

3:0 

1•10"“ 

3.92 • 10”^ 

3 

3.8955 

3•10"® 

0.12 

1 

3.8945 

5•10"® 

0.20 

0.6 

3.8895 

o 

1 

w 

0.27 

0.43 

3.8855 

1■10"^ 

0.39 

0.3 

3.8805 


Table 10. Values of pressure and density ratios and the corresponding position of the zero for a combination of two thermal 
populations with = 8.5 and = 1 keV for different values of ne, 2 . A value rie.i = 3 ■ is adopted here. 


Tie,2 (cm "*) 

P 2 /P 1 

Ue,l/ne,2 

Xo 

1 ■ 10"^ 

1.96 • 10"^ 

3 

3.9005 

1 

0 

CO 

5.88 • 10"^ 

1 

3.8965 

1 

0 

10 

9.80 ■ 10"^ 

0.6 

3.8935 

7 • 10"® 

0.14 

0.43 

3.8895 

1 • 10"^ 

0.20 

0.3 

3.8855 


Table 11. Values of pressure and density ratio and the corresponding position of the zero for a combination of two thermal 
populations with fcsTe,i = 8.5 and fcsTe ,2 = 0.5 keV for different values of ne, 2 . A value ne,i = 3 • is adopted here. 
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Fig. 20. The behaviour of the zero of the total SZ effect produced by a combination of a thermal and non-thermal population 
given in Eq.(^) with ne,rei{pi = 100) = 3 ■ 10“® cm~® as a function of pi (left) and of the pressure ratio Prei/Pth (right). 



n,(pi = 100 ) X 10 ^^ (cm^^) 


Fig. 21. The behaviour of the zero of the total SZ effect produced by a combination of a thermal and non-thermal population 
given in Eq.(^) with pi = 0.5 is shown as a function of ne,rei{pi = 100). 


Pi{s) (see Fig.|22|) around s = 0. At the same time, the total optical depth increases and, as a consequence, the total 
spectral distortion shows a deeper minimum and a higher maximum with respect to the single thermal population case 
(see Fig.lH). It is evident that high values of the density of the cooler population can introduce substantial changes 
to the total SZ effect. 

We also show in Fig.^ the variation of the zero of the total SZ effect as a function of the density of the cooler 
population and of the pressure ratio P 2 /P 1 , respectively. As the cooler population becomes more relevant (i.e., for 
increasing values of its pressure P 2 ), the position of xq moves towards lower frequency values, as expected. In FigJ^ 
the behaviour of xo is shown as a function of the temperature of the cooler population and of the pressure ratio P 2 /Pi 
for a fiducial value of the density He ,2 = 1 • 10“^ cm“^. 
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Fig. 22. The function Pi{s) computed for a combination of two thermal populations with ksTe,! = 8.5 and kBTs ,2 ~ 0.5 keV 
(solid line) and with ksT^,! = 8.5 and kBTe ,2 ~ 1 keV (dashed line) is compared with that of the single thermal population 
with ksTe = 8.5 keV (dotted line). 


6. The spatial distribution of the total SZ effect 

In this section we want to study the effect of the superposition of two distinct electron populations on the spatial 
distribution of the total SZ effect in galaxy clusters. In fact, the thermal and the non-thermal populations have, in 
general, a different spatial distribution as indicated by the spatial extension of the X-ray emission and of the radio 
halo emission in many clusters (see Giovannini et al. 2001). Also, there are indications that the EUV emission in 
some nearby clusters is more extended than the thermal X-ray emission suggesting that there is a warm gas and/or 
a population of relativistic electrons which are spatially more extended than the X-ray emitting IC gas (Lieu et al. 
1999). 

Let us first consider a galaxy cluster in which two different electron populations are present: a thermal population, 
described by a relativistic Maxwellian velocity distribution, and a non-thermal population, with an energy spectrum 
described by a double power-law (see Eq.34). We consider these two populations as independent and spatially super¬ 
posed, consistently with our analysis of Sect. 4. These assumptions are indeed reasonable since in the cluster A2163, 
for instance, the radio halo emitting electrons occupy approximately the same volume of the X-ray emitting, thermal 
electrons (see Eeretti et al. 2001) and, moreover, the relativistic electrons with Eg ~ a few GeV (which produce the 
radio-halo synchrotron emission) do not sensitively affect the IG thermal electrons. 

As for the spatial distribution of the thermal electrons we use here a isothermal 2-iiiodel (see, e.g., Cavaliere & 
Fusco-Eemiano, 1976), according to which the temperature Tg of the IGM is constant and the electron density is 
described by a spherical distribution 

/ ^2 \ - 3 / 3/2 

ne(r) = ngo (^1 -I- , (69) 

where Ueo is the central IG gas density, rg is the core radius of the cluster and (3 = ^mpV^/ksTe is observed in the 
range 0.6 — 1 (see Sarazin 1988 for a review). Under this assumption the angular dependence of the optical depth of 
the thermal population is given by 


t{0) = cjt J ne{r)d£ = arrieo j 


T 

1 + ( - 

Tc 


q -3/3/2 


M 


= neo(JTrcy{9), 


(70) 




Sergio Colafrancesco et al.: Non-Thermal SZ effect in clusters 


29 


0.0010 

0.0005 

0.0000 



Fig. 23. The spectral distortion in units of 2{kBTo)^/(hc)^ produced by a single thermal population with ksT^,! = 8.5 keV and 
rie.i = 3 • (solid line) is compared with that produced by a combination of two thermal populations with kBTe ,2 ~ 1 

keV and ne ,2 = 10“® (dashed line), 3 • 10“® (dotted), 5 • 10“® (dot-dashed), 7 • 10“® (long dashes) and 10“^ (dot-dot-dashes) 


where 9 is the angle between the centre of cluster and the direction of observation. All the angular dependence is 
contained in the function 


Y{e) = 


n¥-\) 

r(f/3) 



(71) 


where 9c = Tc/Da, in terms of the angular diameter distance Da (see Weinberg 1972). 

In the following we will describe in details, for the sake of clarity, the spatial distribution of the approximated 
(to first and second order in r) expression of the SZ effect and we will give also the general expression for the exact 
calculations. 

The angular dependence of the SZ-induced spectral distortion evaluated at first order in r writes as 


Ai{x,9) = T{9)[ji{x) - jo(a;)] , (72) 

where Ai is in units of 2{kBTo)^/(hc)'^ and t{9) is given by Eq.(|70|). Since the spectral distortion can be written in 
the general form Ai{x) = yg{x), the Comptonization parameter writes as 


y{^) = 


(7t 

meC^ 

(7t 

meC^ 

{keTc, 

TTleC^ 
CTj' 


PM = 


(ksTe) / ne{r)M = 


-{9) = 


nieC^ 


PorcY{9) 


(73) 
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Fig. 24. The behaviour of the zero of the total SZ effect for a combination of two thermal populations with fcsTe,! = 8.5 keV 
and kBTe ,2 = 1 keV (solid line) and ksTe,! = 8.5 keV and kBTe ,2 = 0.5 keV (dashed line) is shown as a function of the density 
ne ,2 (left panel) and of the pressure ratio P 2 /A (right panel). A value ne,i = 3 ■ 10~^cm~^ is adopted here. 




kgTg (keV) P2/P1 


Fig. 25. The behaviour of the zero of the total SZ effect for a combination of two thermal populations with fcsTe,! = 8.5 keV 
and ne ,2 = 1 ■ 10“® cm“® is shown as a function of the temperature A:sre_ 2 (left panel) and of the pressure ratio P 2 /P 1 (right 
panel). A value ne,i = 3 • is adopted here. 


where the quantity (ksTg) is defined in Eq.(p^ and depends only on the parameters of the electron distribution (here 
the central pressure is Pq = neoiksTe))- The spectral shape of the SZ effect at an angular distance 6 from the cluster 
center can be written, at first order in r, as 


9{x) 


Ai{x, 9) 

y{^) 


arneoTcY (9) 

7^PorcY{9) 


[ji{x)-jo{x)] 




(74) 
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and it does not contain the angular dependence of the effect. However, at second order in r the spectral distortion 
writes as 


Ai{x,e) = ^T^{d)[j 2 {x) - 2ji{x) + joix)], 


(75) 


and, as a consequence, the second order approximation of the function g(x) is 


g{x,e) = 


1 alniy^Y^e) 

2 ^P,r.Y{9) 

1 (JTmeC^nlQrcY{6) 

2 Pn 


[j 2 {x) - 2ji{x) + joix)] = 


[j 2 {x) - 2ji{x) + joix)] , 


(76) 


which depends explicitly on 9. It is important to notice that the angular dependence of the spectral distortion is 
present both in the Comptonization parameter y and in the spectral shape gix). The dependence of the SZ effect from 
9 is present at all orders n > 1 of the approximations and hence in the full, exact expression of the SZ effect. Thus, 
also for the study of the spatial distribution of the total SZ effect, considering the first order approximation is neither 
appropriate nor consistent. 

For a combination of two electron populations, the total spectral distortion is given, at first order approximation 
in T, by the sum of the single spectral distortions of each population 


Aiixy 9) — ^iihixj 9) T AijiQji^fiix^ 9). 


(77) 


This expression can be written more explicitly as 


Ai(x, 9) 


ythi^)9thix) “t“ ynonthi^)9nonthix) — 
(7 'J' 

- yPo,thrc,xYthi9)gthix) + 

rrieC 
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(jj' 


' Po,nonthT"c,rad^nonth (^) — 


rrieC^ 

(T 'T~' 

- :^Po,th[rc,xYthi9)gthix) H 

rrieC^ 

F Xc.radPYnonth i9)gnonth (^)] 


(78) 


where P = Po,nonth/Po,th is the pressure ratio at the cluster center. We assume here that the spatial distribution of 
the thermal and of the non-thermal electron populations are given by a /3-model with core radii Tc^x and Xc^rad and 
parameters f3x and I3rad, respectively. Such an expression is also valid in the case of the combination of two thermal 
populations. At higher approximation orders in r, the total spectral distortion contains cross-correlation terms of the 
parameters of each population and thus it is not possible to write it as a linear combination of the spectral distortions 
of the separate electron populations. An exact derivation (or at least the 3’’'^ order approximation in r) is required 
to describe correctly the spatial dependence of the total SZ effect in galaxy clusters. The general calculation of the 
spatial distribution of the total SZ effect can be made using the exact expressions derived in Eqs.(21-23, 62 and 63) 
in which the spatial dependence of the thermal and non-thermal electron populations are considered. 

In Figj2^ we show the spatial behaviour of the total SZ effect in the case of a Coma-like cluster. We use specifically 
the following parameters: Xf- x = 0.42Mpc, Tetrad = 0-4Mpc, Rx = 4.2Mpe, Rhaio = 1-25Mpe, 
I3x = 0.75, Prad = 0.8, neo,x = 2.89 • 10“^ /i|q cm“^, k^Tf. = 8.21 keV. The non-thermal electron population fitting 
the radio halo spectrum is taken with spectral parameters pi = 1000 and a = 2.5. 

The total SZ effect shows a peculiar spatial behaviour which is particularly evident at the frequencies x ^ 2.3 where it 
attains its minimum value and at frequencies x ^ 6.5 close to its maximum. At a: ^ 2.3 the spatial profile of the total 
SZ effect is predicted to be higher than that due only to the thermal population because the non-thermal SZ effect 
sums up in the region of the radio halo which extends up to ^ l.Sh'^o Mpe (see Fig.^^. At larger distances from the 
cluster center, the SZ effect is dominated by the thermal component and resemble its spatial behaviour. An opposite 
behaviour is expected at frequencies higher than that of the zero of the total SZ effect. For instance at a; ~ 6.5 the 
spatial profile of the total SZ effect is lower than that of the single thermal effect because the spectrum of the total SZ 
effect is lower than that of the thermal case. So, a depression in the spatial profile of the total SZ effect is expected at 
radii smaller than the radio halo extension Rhaio- At distances larger than the radio halo extension, the total SZ effect 
is again dominated by the thermal component. The amplitude of the non-thermal SZ contribution depends from the 
pressure ratio P = Prei/Pth as shown in Fig.^. Such spatial features of the total SZ effect should be detectable with 
sensitive instruments with good spatial resolution and narrow spectral frequency band. Such instrumental capabilities 
can also be able to detect the transition between the ambient thermal SZ effect and the non-thermal SZ effect expected 
in radio relics and in the radio lobes of active galaxies found in the cluster environment. 
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Fig. 26. The spatial dependence of the total SZ effect produced by the combination of the thermal and relativistic electrons in 
a Coma-like cluster. We show the thermal SZ effect (dots) and the total SZ effect with P = 0.05 (solid), 0.49 (short dashes) and 
1.48 (long dashes). Calculations are done for an a-dimensional frequency of a: = 2.3 where the total SZ effect has its minimum 
value (upper panel) and for x = 6.5 where the total SZ effect has its maximum value (lower panel). 

7. Application to specific galaxy clusters 

The non-thermal SZ effect may already have been observed as a part of the SZ signal from some hot galaxy clusters 
which possess also non thermal phenomena like radio halos or relics. 

In the Coma cluster, in fact, Herbig et al. (1995) detected a strong SZ effect at the level of y ~ lO”"* from an 
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observation at low frequencies Vr ^ 30 MHz. Since the Coma cluster has a bright and extended radio halo (see, e.g., 
Deiss et al. 1997), it is well possible that some fraction of the SZ signal might be contributed by the non-thermal 
SZ effect. However, to disentangle the non-thermal SZ effect from the thermal one a detailed spectral coverage with 
reasonably good sensitivity in the frequency range from x ~ 2.5 up to x ^ 8 is needed, with particular care to the 
region x ~ 3.8 — 4. At the time of the acceptance of this paper, the available SZ observations for Coma at these 
relevant frequencies are still lacking and we will discuss this case in more details in a further paper (Colafrancesco & 
Marchegiani 2002). 

In the following, nonetheless, we will discuss specifically the case of A2I63 for which there are observations at different 
frequencies which cover the interesting part of the SZ spectrum and can be used to put constraints on the presence 
and on the nature of the non-thermal SZ effect in this radio-halo cluster. 

7.1. The case of A2163 

The cluster A2163 {z = 0.203) is one of the hottest {ksTe ^ 12.5 keV) clusters which possesses a giant radio halo 
(Herbig & Birkinshaw 1994; Feretti et al. 2001) with diameter size of ^ (2.9 ± 0.1 )/i^q^ Mpc which is centered on the 
peak of the X-ray emission. The slope of the power-law synchrotron spectrum is ~ 1.6 ±0.3 and has been estimated 
from radio data at 1.365 and 1.465 GHz (Feretti et al. 2001); this is consistent with the results obtained by Herbig & 
Birkinshaw (1994) in the frequency range 10 MHz - 10 GHz. There is, however, no evidence of hard X-ray excess due 
to a non-thermal component in the BeppoSAX data of this cluster (Feretti et al. 2001). 

A strong SZ effect has been observed in A2163 at different frequencies. In particular the SZ effect spectrum has 
been observed from BIMA at 28.5 GHz (LaRoque et al. 2002), from DIABOLO at 140 GHz (Desert et al. 1998) and 
from SuZIE at 140, 218 and 270 GHz (Holzapfel et al. 1997; these data are dust-corrected in LaRoque et al. 2002) thus 
including and bracketing the null of the thermal SZ effect. Gonsidering a pure thermal SZ effect, the previous data on 
A2163 are fitted with a Gompton parameter yth = (3.65 ± 0.40) • 10“^ and with the addition of a kinematic SZ effect 
whose amplitude corresponds to a positive peculiar velocity Vp = dlhl^gg km s“^ (Carlstrom et al. 2002) whose large 
uncertainties, however, makes it consistent with a zero value. Note that, in principle, the addition of a positive-velocity 
kinematic SZ effect contributes to provide a deeper negative SZ signal in the region between the minimum and the 
zero of the SZ effect. 

Here, we re-analyzed the data on A2163 trying to put constraints on the possible presence of a non-thermal SZ 
effect by fitting the available data with a combination of a thermal and non-thermal SZ effect. 

The thermal population in A2163 has a temperature of ksTe = 12.4±0.5 keV and a central density of ne,th — 6.82- 
10“^ cm“^. The parameters describing the spatial distribution of the IC gas are Xc = 0.36 h'^Q Mpc and /3 = 0.66 (Elbaz 
et al. 1995; Markevitch et al. 1996). With these values the optical depth towards the cluster center is TQ^th = 1-56 • 10“^ 
and the central Gomptonization parameter is = 3.80 • 10“^. 

We show in Fig. the relativistic, thermal SZ effect expected for A2163 compared to the available data. The 
thermal SZ effect fits the data with a = 1.71 and hence it is statistically acceptable. However, we will show in the 
following that the inclusion of a non-thermal component of the SZ effect is able to improve sensitively the fit to the 
data. 

We consider a non-thermal population with a double power-law spectrum with parameters a\ = 0.5, oti = 2.5, 
•pcr = 400, p 2 ^ oo and we set pi and the density ne,rei{pi), with pi = 100, as a free parameters. We also assume 
that the spatial distribution of the non-thermal population is similar to that of the thermal population as indicated 
by the extension of the radio halo in A2163 (see, e.g., Feretti et al. 2001). In our calculations, we fix pi and we search 
for the value of rie^reiiPi) which minimizes the We show in Table ^ and in Fig. n the results of our calculations. 
We see that the data are better fitted with values of pi which correspond to the electrons with pi ^ 10^; we see also 
that these electrons can produce a detectable non-thermal SZ effect with a value of the pressure {Prei/Pth ^ 0.3) and 
the density ratio n ^ 270 which do not imply a strong influence on the dynamical and thermal state of the IC gas. 
A model in which pi = 100 and ne,rei{pi) ~ 2.5 • 10“® cm“^ best fits the data with a Xmm = 1-0534, much lower 
than the Xmin — l-l’l yielded by the single thermal population. The previous parameters point to a non-thermal 
population which carries a non negligible pressure contribution P^ei ~ 0.29Pth and which has a spectrum not extended 
at momenta lower than Pmin — 100, corresponding to Emin ~ 50 MeV. 

This result does not imply that the fit with a combination of thermal plus non-thermal populations is statistically 
excluding the single thermal population model. However, our analysis shows that: i) a non-thermal SZ effect is produced 
by relativistic electrons producing radio halo emission has to be present in A2163. The analysis indicates that its 
amplitude could be appreciable (see Fig.27) and corresponds to a pressure in relativistic particles Prei ~ O.SPth', H) 
the possibility of having SZ observations with better precision can offer the possibility to disentangle between the 
thermal and any non-thermal component of the SZ effect; in) the detection of a non-thermal SZ effect can set strong 
constraints on the nature of the non-thermal population and on its feedback on the thermal one. 
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In these respects, it is appealing that the physical characteristics of the non-thermal population can be constrained 
through a detailed study of the SZ effect observed in the same galaxy cluster. Once the value of pi has been set, 
the quantity Ue^rei {Pi ) can be constrained with a quite good precision and vice-versa. This example shows one of the 
potential uses of the SZ effect to obtain information on the properties of different electronic populations which are 
residing in the atmospheres of galaxy clusters. 


Pi 

Tie,relipl) (cm "‘) 
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0.1 
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1.0536718 
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10 “® 

0.25 
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0.26 
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1.0534818 

100 
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0.29 
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1.0534685 

1000 

7.991 • 

10 ~^ 

0.38 

270 

1.0534679 

10000 

2.527 • 

10~2 

0.46 

270 

1.0534679 


Table 12. For each value of pi are indicated the value of the density ne,rei{pi ~ 100) which provides the best fit and the 
corresponding values of pressure (P = Prei/Pth) and density (n = ne,th/ne,rei) ratio with the corresponding value of for the 
total SZ effect in A2163 produced by a combination of a thermal population and a non-thermal double power-law population. 
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Fig. 27. Theoretical expectations for the spectrum of the SZ effect in A2163. We show the fit to the available data yielded by 
a thermal population (solid curve) and the expectations obtained from a combination of thermal and non-thermal populations 
with Pi = 100 for a value of the pressure ratio Prei/Pth = 0.29 (dashed curve), which provides the best fit. 


Using our general approach described in the previous sections, we also considered the calculation of the total SZ 
effect produced by a combination of two thermal electron populations, with the warm population temperature in the 
range ksT ^ 0.1 — 1 keV and density in the range rig ,2 ^ 10“^ — cm“^. We also considered two different spatial 

distributions of the hot and warm populations: in the first case the two populations have the same spatial extent 
and in the second one the cooler population is spatially more extended than the hotter one. In this last case we use 
parameters rc ,2 = l-5rc,i and (32 — 0.5. The results of the fit to the data of A2163 are reported in Tablep^ and in 
Fig.^, and in Table|l^ and Fig.^. The available SZ observations allow also in this case to set constraints on the 
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fcsT’e .2 (keV) 

Tie,2 (cm "*) 

Pi/Pi 

^e, 1 /^e,2 


0.1 

3.922 • 10“^ 

4.64 • lO""" 

0.17 

1.5633934 

0.5 

8.508 • 10"® 

5.03 • 10~2 

0.80 

1.2395049 

1 

4.222 • 10"® 

4.99 • 10"^ 

1.61 

1.2307970 


Table 13. For each value of fcsTe ,2 are indicated the value of the density n 2 which provides the best fit and the corresponding 
values of pressure and density ratio with the corresponding value of for the total SZ effect in A2163 produced by a combination 
of a two thermal populations with the same spatial distribution. 


parameters ne ,2 and Te ^2 A best fit which minimizes the Xmin i® obtained for values fc_BTe _2 = 0.5 — 1 keV, while the 
case kBTe ,2 = 0.1 keV seems to be disfavoured by the data. 

Finally, we want to remark that, in the set of models we considered in this section, the best fit to the data (with 
a Xmin ~ 1.05) is obtained from a combination of thermal and non-thermal populations, a case which seems to be 
favoured by the present data on A2163. Future SZ observations of A2163 with higher spectral resolution and sensitivity 



X 


Fig. 28. The spectrum of the SZ effect in A2163 obtained from a combination of two thermal populations with fesTe,! = 12.4 
keV and kBTe ,2 ~ 0.5 keV and with the same spatial distribution. We show the cases of P 2 /Pi'- 0 (solid line), 1.77• 10“^ (dashed 
line), 5.03 ■ 10“^ (dotted line, which yields the best fit in this case), 5.91 • 10~^ (dot-dashed line). 


fcsTe ,2 (keV) 

ne ,2 (cm "*) 

Pi/Pi 

ne,l/ne,2 


0.1 

1.624 • 10~^ 

1.92 • lO”"" 

0.42 

1.5911710 

0.5 

3.531 • 10"® 

2.09 • 10"^ 

1.93 

1.2504161 

1 

1.752 • 10"® 

2.07 • 10"^ 

3.89 

1.2412285 


Table 14. For each value of kBTe ,2 are indicated the value of the density n 2 which provides the best fit and the corresponding 
values of pressure and density ratio with the corresponding value of for the total SZ effect in A2163 produced by a combination 
of a two thermal populations with different spatial distributions. 
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Fig. 29. The spectrum of the SZ effect in A2163 obtained from a combination of two thermal populations with fcsTe,! = 12.4 
keV and fcsTe ,2 = 0.5 keV and with different spatial distributions as explained in the text. We show the cases of P 2 /Pi- 0 (solid 
line), 5.91 ■ 10“® (dashed line), 2.09 • 10'^ (dotted line, which yields the best fit), 2.96 ■ 10~^ (dot-dashed line). 


will allow to set stronger constraints on the non-thermal population. In particular, the possibility to measure accurately 
the frequency location of xo could yield direct information on the value of the pressure in relativistic particles, Prei, 
confined in the cluster atmosphere and on their energy distribution. 

8. Discussion and conclusions 

Since many galaxy clusters with available SZ measurements show also the presence of an extended radio halo and 
other non-thermal phenomena, it is relevant to understand quantitatively the relevance of the non-thermal SZ effect 
in galaxy clusters and to disentangle it from the true thermal SZ effect which is extensively used for cosmological 
studies. It is also relevant to use the limits on the non-thermal SZ effect to put constraints on the spectrum of the 
non-thermal population which is responsible for the high-energy phenomena occurring in galaxy clusters. This proves 
to be a powerful and complementary analysis of the specific studies of non-thermal phenomena carried out at other 
wavelenghts. 

In such a context, we have shown in this paper that the spectra of the thermal and non-thermal SZ effects are 
distinctly different and that the overall spectrum of the SZ effects measures the energy densities in the thermal and in 
the non-thermal electron populations, separately. We have calculated here the spectral shape of the non-thermal SZ 
effect in the Thomson limit using an exact approach which makes use of the full relativistic formalism and goes beyond 
the Kompaneets and the single scattering approximation. Moreover we evaluated, for the first time in a consistent way, 
the total SZ effect arising from a combination of different electron populations residing in the same galaxy cluster: 
we considered specifically the cases oi i) a thermal plus a non-thermal electron populations and ii) two thermal 
populations with different temperatures and densities. Such a derivation can be obtained, nonetheless, for an arbitrary 
set of different electron populations. 

We have shown in our approach that the spectral distortion of the CMB due to both the thermal and the non- 
thermal SZ effect can be written in a general, single form which depends separately on the Comptonization parameter 
and on the spectral shape of the specific SZ effect. Such an expression is valid for both thermal and non-thermal SZ 
effects produced by single electron populations as well as for the overall SZ effect produced by a combination of different 
electron populations. We have also shown that the spectral dependence of the SZ effect, contained in the function g{x), 
cannot be separated from the dependence on the cluster parameters even at the first order approximation in r, where 
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g{x) depends on the cluster (effective) temperature. At higher-order approximation in r, g{x) depends, moreover, also 
on T. In the exact approach, the spectral shape of the SZ effect depends strongly on the cluster parameters themselves 
and specifically from the pressure and the optical depth of the considered electron population. 

We have calculated in details the precision at which each approximation order in t reproduces the exact form of 
the SZ effect: we have shown that the third order approximation in t yields a precision ^ 1% at each interesting 
frequency for both the non-thermal and the thermal SZ effect. 

We have shown that the main spectral differences for the non-thermal SZ effect with respect to the thermal SZ effect 
are a maximum of the spectral shape which is moved to higher and higher frequencies for increasing values of the lower 
cutoff in the momentum distribution, pi. As a consequence, also the null of the non-thermal SZ effect, xq, is moved to 
higher frequencies for increasing values of pi. Since the pressure Prei of the non-thermal electron distribution depends 
on Pi , the position of xq decreases with increasing values of Prei for a fixed normalization of the non-thermal electron 
density. 

We have shown in this paper that the shape of the exact non-thermal SZ effect is different from the 1®* order 
approximation estimate used so far by different authors (Birkinshaw 1999, Ensslin & Kaiser 2000, Blasi et al. 2000, 
Shimon & Rephaeli 2002). We have shown that the 3’’'^ order approximation is required to reproduce the exact results 
within a ~ 1% precision at all relevant values of x. 

We have also shown that the SZ effect from a single non-thermal population in clusters with non-thermal phenomena 
is not a realistic description of the effect. So, we have derived for the first time a consistent expression for the total 
SZ effect produced by a combination of a thermal and a non-thermal populations, like is the realistic case in many 
galaxy clusters. The shape and the amplitude of the total SZ effect is quite different from that suggested in previous 
papers (Ensslin & Kaiser 2000, Blasi et al. 2000). The main difference in the derivation of the overall SZ effect is 
that it is not just given by the sum of the separate thermal and non-thermal contributions (as is the case for the first 
order approximation in r) but it is given by a properly weighted (according to the optical depths of each separate 
population) combination of the thermal and non-thermal contributions, as discussed in Sect. 4. Also, we considered in 
our approach the effect of multiple scattering and higher (n > 1) order approximation terms in the analytical derivation 
of the overall SZ effect. Moreover, we also discussed not only the case of i) a single power-law electron spectrum (as 
was done also by the previous authors) but also the more consistent case of a ii) double power-law electron spectrum 
which is able to fit both the radio halo and the hard X-ray spectra in clusters. We showed that, while in the case 
i) low-energy electrons are the dominant source of non-thermal Compton scattering and produce an increase of SZ 
signal in the region x ^ 7 (where the thermal SZ effect has its maximum), in the case ii) high-energy electrons are 
the dominant source of non-thermal SZ effect and produce very large shifts of CMB photon frequencies producing a 
deprivation of photons also in the region of the maximum of the thermal SZ effect (x ~ 7). 

A relevant result of our analysis shows that the location of the zero of the total SZ effect depends on the pressure 
ratio P = Prei/Pth between the relativistic and thermal pressures of the two different electron populations. It increases 
non-linearly with P up to values xq ~ 4.2 for P up to values « 1.5. This is a unique and relevant feature of the non- 
thermal SZ effect in clusters since it yields, in principle, a direct measure of the total pressure in relativistic non-thermal 
particles in the cluster atmosphere, an information which is not easily accessible from the study of other non-thermal 
phenomena like radio halos and/or EUV or hard X-ray emission excesses. Thus, the detailed observations of the 
frequency shift of xq in galaxy clusters provides unambiguously a constraint to the relativistic particle content in the 
cluster atmosphere. In addition, such a measurement is crucial to determine the true amount of kinematic SZ effect in 
clusters since it is usually estimated from the residual SZ signal at the location of the zero of the thermal (relativistic) 
SZ effect. Due to the steepness of the SZ spectral shape in the region of the null of the overall SZ effect, the possible 
additional non-thermal SZ signal must be determined precisely in order to derive reliable limits for the kinematic SZ 
effect. This fact can severely limit the possibility to measure effectively the kinematic SZ effect in radio-halo clusters. 
An accurate subtraction of the non-thermal SZ effect needs a specific observational strategy with measurements in at 
least three frequency ranges x~2 — 3, 3.8 — 4.5 and ~ 6 — 8. 

We verified that our approach is consistent with the exact derivation of the monochromatic redistribution function 
containing both the relativistic and quantum corrections already given by Fargion et al. (1997) and Sazonov & Sunyaev 
(2000) in the Thomson limit and in the single scattering limit. In the Thomson limit, nonetheless, we have generalized 
the full derivation of the SZ effect to the case of multiple scattering and to any approximation order in r. 

After the submission of our paper, Shimon & Rephaeli (2002) presented a computation of the SZ effect produced 
by non-thermal electron populations in a few clusters like Coma and A2199. However, we want to stress that their 
derivation of the non-thermal SZ effect is still approximate since they consider only the first order approximation in 
T (in which case the total SZ effect is the sum of the thermal and non-thermal effects) and the single scattering limit. 
We have shown in our paper that such an approximated description is neither complete nor precise since at least the 
third order approximation calculation in t is required to reproduce the exact derivation of the total SZ effect within 
~ % accuracy level. 
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We have also generalized our derivation of the total SZ effect in galaxy clusters to the case of a combination of 
different thermal electron populations (see Sect. 5). Any additional cool IC gas component produces a tightening of the 
photon redistribution function, an increase in the total optical depth and hence a substantial change in the spectral 
distortion at the minimum and at the maximum of the SZ effect. The location of the zero of the SZ effect xq also 
decreases in frequency due to the presence of the cooler component which decreases the pressure ratio P 2 /Pi- Thus, 
the possible detection of an additional cold component in the cluster atmosphere through observations of the total 
(thermal plus thermal) SZ effect allows to test the possible thermal origin of the EUV excess observed in several nearby 
clusters (Lieu et al. 1999, 2000). Our analysis of the cluster A2163 does not indicate a relevant role of the warmer 
component to the total pressure of the cluster atmosphere, consistently with the present indication of the absence of 
a strong EUV excess in this cluster. We will address the specific analysis of other clusters elsewhere. 

Beyond the relevance of the study of the non-thermal SZ effect as a bias for the cosmologically relevant thermal 
SZ effect, the non-thermal SZ effect has also a crucial astrophysical relevance as a barometer for the presence of any 
additional population of electrons with both a non-thermal or a thermal energy spectrum. In fact, the non-thermal 
SZ effect actually measures the total pressure of the non-thermal electron population and hence yields constraints to 
its energy spectrum. Analogously, the additional thermal SZ effect produced by a cooler thermal component provides 
information on its temperature and density. Specifically, we have shown that combining information from the cluster 
radio-halo observation and non-thermal SZ measurements, one is able to determine the shape and the extension of the 
relativistic electron spectrum. Using the available observation of one of the best example of radio-halo clusters with 
multifrequency observations of the SZ effect, A2163, we applied our method and derived a limit on the low-energy 
cutoff of the relativistic electron spectrum of Emin ~ 50 MeV, as well as constraints on the momentum spectrum 
which is required to be flat enough (/e ~ p“°-^) below p ~ 400 to avoid destructive feedback effects on the thermal IC 
gas. 

Relativistic electrons with such flat energy spectrum do not produce a relevant extra heating and/or extra X-ray 
emission with respect to the thermal IC gas in A2163. Also the IC energy losses of such relativistic electrons do not 
yield a substantial hard X-ray emission, in agreement with the available limit on A2163 obtained from BeppoSAX 
observations. 

We have finally shown that the presence of a non-thermal SZ effect also influence the spatial profile of the total SZ 
effect of a typical radio-halo cluster. In fact, the region occupied by the radio halo shows an increment (decrement) 
of the signal at frequencies near the minimum (maximum) of the SZ effect. This happens as a consequence of the 
difference in the spectral shapes of the non-thermal SZ effect with respect to the thermal one. 

The specific spectral and spatial features of the non-thermal SZ effect makes it possible to detect it through a 
multifrequency observation with high sensitivity and narrow-band detectors. The optimal observational strategy is to 
observe galaxy clusters in the frequency range a: ^ 2 — 8 where the peculiar spectral features allow clearly to disentangle 
the non-thermal SZ effect from the thermal one. The PLANCK surveyor experiment has the capabilities to detect and 
map the non-thermal SZ effect in a large number of nearby radio-halo clusters. However, dedicated experiment with 
high sensitivity and narrow band spectral coverage are also adequate to detect the non-thermal SZ effect in radio-halo 
galaxy clusters. 


Acknowledgments. We thank the Referee, N. Itoh, for several useful comments and suggestions which allow us to 
improve the presentation of our results. 

Appendix A: The total SZ effect for a combination of two electron populations: analytic 
approximations 

In this Appendix, we derive detailed expressions for Itot{x) which are approximated at first and second order in r 
according to the general expression: 

I{x) = Iq{x) + t[Ji{x) - Joix)] + ]^t'^[J2{x) - 2Ji{x) + Joix)]. (A.l) 

We remind the reader that the function Ji{x) is given by Eq.(18): 

/ + 00 

Io{xe~‘^)Pi{s)ds, (A.2) 

-OO 

where Pi{s) is given by Eq.(^): 
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Inserting Eg. (p2[) in the previous Eg. for Pi(s), we obtain: 

poo poo 

Pi{s) = CA dpfA{p)Ps{s;p)+CB / dpfB{p)Ps{s;p) 
Jo Jo 

= c^-Pia(s) + cbPi_b(s). 


(A.4) 


Inserting Eg.(A^) in Eg.(|A.2|), we also obtain the expression: 


/ -hoo p-\-co 

Io{xe~'')PiA{s)ds + cb / Io{xe~'')PiB{s)ds 

-oo j —oo 

= caJia{x) + cbJib{x). 

Inserting egs.(A.5), d^) and (p4[) in Eg.(A.l), we derive, at first order in r, the expression: 


(A.5) 


= ('^A + TB)[cAJlAix) + CbJib{x) - Mx)] 

= taJia{x) + tbJib{x) - TA-h{x) - TbJo{x) 
= ta[Jia{x) - Jo(a:)] +tb[Jib{x) - Jo(a;)]. 


(A.6) 


As expected, the SZ effect at first order in r is given by the sum of the separate SZ effects produced by the two electron 
distributions, respectively. 

To evaluate the effect at second order in r, it is necessary to calculate the expression: 


/ + 00 

/o(xe“'*)P2(s)ds. 

-OO 

The function P 2 {s) can be derived by Eg. (i using Eg. (|A.4D : 

P 2 {s) = Pi(s)®Pi(s) 

= [caPia{s) + cbPib{s)] ® [caPia{s) + cbPib{s)\ 

= c\Pia{s) 0 Pia{s) + 2cacbPia{s) 0 Pib{s) + c%Pib{s) 0 Pib{s) 
= c\P2Aa{s) + 2caCbP2Ab{s) + c\P2Bb{s)- 


(A.7) 


(A.8) 


Inserting this expression in Eg.(A.7), we obtain: 

/ + 00 |■ + CO 

Io{xe~'')P2AA{s)ds+ 2 cacb / Io{xe~‘^)P2AB{s)ds + 

-oo J —oo 


/ + 00 

Io{xe~‘^)P 2 BB{s)ds 

-oo 

= c\J2Aa{x) + 2caCbJ2Ab{x) + C%J2Bb{x). 

From Eg. B) we derive the expression of the second order correction to the total distorted spectrum 
~ 2^'’'^ Xb)‘^\<?A'^2,AA + ‘2caCbJ2AB + C%J2BB ~ ‘2-CaJiA + 

—2cbJib + Jq] = 

= 2^x\J2aa + ‘^taTbJ2AB — ‘^xa{ta + tb)Jia — 2tb(t^ + tb)Jib + 

+ (t^ + 2tatb + = 

= ■;::{xA[d2AA — 2JiA + > 70 ] + Tb[J2BB ~ ‘2-JlB + Tq] + 


(A.9) 


+2tatb[J2Ab — JiA — Jib + -^o]}. 


(A.IO) 


Notice that at second order in t there is an additional term describing the probability that a CMB photon can suffer 
first a scattering from the electrons of the distribution /a and then another scattering from the electrons of the 
distribution /b- 

In a similar approach, we can derive the further corrections to the total distorted spectrum evaluated at any order 
in r; these will contain a larger and larger number of cross-scattering terms. In a completely analogous way, we can 
also consider the case of three and more electron distributions. 
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